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FOREWORD 


The Industrial Mathematics Society is a professional organiza- 
tion whose objective it is to extend the understanding and application 
of mathematics in industry. Founded early in 1949, it marked the 
first organized attempt of any group to foster a closer relationship 
between mathematics and industry. 

It promotes the cause of mathematics in industry through month- 
ly meetings at which technical papers are presented, through 
lectures by nationally prominent scientists and engineers, through 
panel groups which concentrate in specialized fields of knowledge, 
and through publication of worthwhile papers in its annual volume. 

Membership is open to engineers, scientists, and mathematicians 
who are concerned with the effective use of mathematics in industry. 
Annual dues are $5.00 per year including a subscription to the 
annual volume, Industrial Mathematics. Application forms and 
information may be obtained by writing to the Secretary, Industrial 
Mathematics Society, 100 Farnsworth, Detroit 2, Michigan. 
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The Engineering Applications of Statistics 


By David N. Collins 
Chrysler Corporation 
Engineering Division 
Detroit 


The operation of any business, whether it is building automobiles 
or building a hamburger, brings an endless succession of questions. 
Can I use a different steel? Can it be lighter and still stand the 
strain? Did the butcher take me on that last batch of meat? The 
businessman’s feeling of utter helplessness which accompanies any 
serious look at the day's list of questions has welcomed the help of 
engineering and research groups, or, lacking them, anyone at all 
who claims to have The Answer—any answer. 

But engineers, to answer the questions of business, have ques- 
tions of their own. How many service failures will occur with this 
new steel? If this part fails at 4,307 pounds, will the other 999, 999 
parts we plan to make also fail at 4,307 pounds, exactly? Are there 
any standards for hamburger? These questions get mathematical, 
and the branch of mathematics which can answer them best is sta- 
tistics. It is the engineer's feeling of helplessness which adds the 
mathematical statistician to the payroll. 

The addition of a statistician to the engineering staff may multi- 
ply the engineer's problems, rather than divide them. "What isa 
regression? Why won't you plot a smooth curve to show our data? 
What if it does have nine inflection points?" 

The statistician has his own problems. "Don't these people 
understand statistical significance? How can they expect conclusive 
results from only one test, or even two?" 

To sum up, there is a gap between the engineer and the statis- 
tician which must be bridged before the statistician is of much value 
to engineering. This paper will attempt to define the gap, and to 
suggest several bridges through which useful results can be ob- 
tained. 


ENGINEERING ENDS AND STATISTICAL MEANS 


Engineers use mathematics. Statisticians use mathematics. But 
the mathematics of the engineer is applied mathematics, tied to ob- 
servations of physical quantities, and it must be resolved into phys- 
ical effects. Statisticians use artificial arrays of numbers which 
may (or may not) be associated with physical observations. Luckily 
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these arrays and other techniques can greatly expand the range and 
sensitivity of engineering experiments. But it does not follow that 
every statistical operation has a physical counterpart. 

ENGINEERS must make accurate and reproducible measurements, 
characterize and compare test methods, forecast responses to ser- 
vice conditions. 

STATISTICIANS can separate and measure variances of groups of 
observations, measure the probability of real differences (not due to 
chance), evaluate constants in an equation of correlation. 

The engineering applications of statistics are the use of the sec- 
ond set of functions for the benefit of the first. 


ENGINEERING RESULTS FROM STATISTICAL METHODS 


Probably the best non-mathematical link between statistics and 
engineering is through a discussion of individual applications for 
each method in turn. One of the simplest is the separation of var- 
iances within groups of observations; this is called the analysis of 
variance. 


THE ANALYSIS OF VARIANCE 


This is probably the most powerful single tool of statistics as 
applied to engineering. To illustrate, if a number of observations 
(of anything), made under varying conditions, are arranged in some 
kind of rectangular array, the result would be the vertical bars of 
Figure 1. These observations have an overall average (the plane of 
Figure 2) but are also divided by the rows and columns of the array 
into smaller sub-groups. each with its own average (Figure 3). 
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Figure 1 Figure 2 Figure 3 


Each of these averages is accompanied by its calculated variance, 
which measures the variability of the observations represented by 
the average. So far, the arrangement of the array has been arbi- 


trary, without relation to any physical quantity. But, if the rows 


and columns can be made to stand for physical conditions, then sta- 
tistics and engineering come together with beneficial results that 
can be obtained in no other way. These quantities can be calculated: 
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1. The effect of each of the independent variables (usually one, 
two, or three) by which the observations were classified into rows, 
columns, and other groups. 

2. The effect of each combination of variables taken two, three, 
or more, at a time. 

3. The residual variation not accounted for by such a classifica- 
tion. 

4. The probability that each of the effects found is due only to 
chance variation of the data. On this probability one can decide 
whether or not these effects are real and be right, say, nineteen 
times in twenty. In a word, the variance of all the observations has 
been analyzed into one sectidn for each independent variable and 
combination of variables, with an unexplained remainder against 
which the others may be judged for chance differences. (Figure 4) 


ROWS 








COLUMNS 








INTERACTION ——— 





RESIDUAL 








TOTAL 36 TOTAL 
Figure 4 


If the rectangular array is completely filled, that is, if there is 
one observation for each intersection of row and column, then this 
is a factorial design. There is an observation for each combination 
of the conditions assigned to the rows and columns. This is the sta- 
tistical design which is most often compared with the classical de- 
sign, which changes one variable ata time, while holding all the 
others constant. A detailed comparison of this type has been made 
by Duffett. The advantage of the factorial over the classical design 
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is that less than twice as much work gives about four times as much 
information. 


TESTING THE REALITY OF A DIFFERENCE 


Probability calculations over the difference between two averages 
do not have the mathematical glamor of an analysis of variance 
(where "one problem makes a publication"), but they do provide a 
mathematical basis for every decision based on the reproducibility 
of resuits. The only numerical data required are the two averages 
and the variances associated with each average. If the difference 
between the averages is expressed in units of the combined variance, 
a figure can be found to express the probability that the difference is 
due only to chance variation, and that further measurement would 
show that the averages are not really different at all. Such a calcu- 
lation is known as the application of student's ''t'', which is the unit 
of measurement in which the difference is calculated. 

Practically, the use of calculated probabilities like these is 
limited to cases in which the difference between averages is not 
obvious by inspection. But the difference by inspection is usually 
not as obvious as it looks. Of the two columns below, the second 
shows an average increase of 61.9 percent over the first. 


4. 48 11.01 
3.59 5. 27 
6. 42 7.16 
Av. 4.83 7.82 Increase 61.9% 


Very straightforward and respectable. But the calculation of a prob- 
ability figure through student's "'t'’ shows that the probability of a 
real difference is only 0.66, that is, that one out of three samples 
of this size would show a difference this large if there were no real 
difference at all. 

Naturally, this is a loaded example. But it is loaded no worse 
than much laboratory data from which differences are inferred and 
broad conclusions drawn. Remember that this “proof” of a differ- 
ence still has a 66 percent chance of being confirmed by more data, 
but there is also a 33 percent chance that no difference exists. How 
much money would you bet on it? 

Two things give this data a false appearance of authority. First, 
the results are calculated to two decimals, in spite of a variance of 
two to four units in each group. Second, the percentage increase 
looks impressive, although three samples which disagree so widely 
are too few to give a percentage difference any validity. Probability 
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calculations distinguish between guesses and proof. 


CORRELATION 


Statistics applied to graphical data produces measures of corre- 
lation. From any group of x, y observations the line or curve of best 
fit can be calculated, as well as a number which is a measure of the 
closeness of fit (Figure 5). The line or curve is calculated by the 


CORRELATION 
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Figure 5 


least squares method, which minimizes the sum of squares of the 
distances from each point to the calculated line. These are mathe- 
matical operations, of course, and the calculated line need not have 
any relation to a physical property of the data. A straight line, ora 
curve of any degree may be arbitrarily imposed on the data. 

In the equation, 


y = ax+b, 


the calculated values are the slope, a, andtheintercept, b. Although 
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the curves can be fitted mathematically, it is unusual to find any- 
thing but the straight line fitted in practice. It is easier to trans- 
form y into y, log y, or some similar function, and to fit the 
straight line than to go through the calculations for second degree 
or higher curves. 

Another useful correlation which can be calculated is 


y = ax; + DX, + CXs3..., 


which gives the relative effects of a number of independent variables 
on a single result (Figure 6). This is most useful for measuring the 
relative importance of several variables on a single physical prop- 
erty, such as the effects of the alloying elements of steel on its ten- 
sile strength. 


MULTIPLE CORRELATION 














Figure 6 


This multiple correlation is more quantitative than the analysis 
of variance which was mentioned earlier; where the usual analysis 
of variance only identifies the significant variables, multiple corre- 
lation evaluates their relative importance. Multiple correlation re- 
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quires a larger mass of data than the analysis of variance, and the 
calculations are more complex; however, the requirements for pre- 
setting the conditions of measurement are not as strict. The anal- 
ysis of variance more easily shows interaction among the variables, 
which is only indicated in multiple correlation by a poor fit at some 
points between the observed data and the calculated line. 


APPLICATION 
DESIGN OF EXPERIMENTS 


Not many technical men, these days, deny that the use of sta- 
tistics improves the quality of research work. Most problems of 
application arise from the fact that the statistician is only called in 
after the laboratory work is done. Without pre-set conditions of 
measurement, statistics cannot be used to best advantage. The pre- 
set conditions are never reached without prior planning. 

In the design of experiments, the statistician acts best as a 
referee, advising and arbitrating between two more important peo- 
ple—the man who will do the laboratory work and the man who asked 
the question which the work is designed to answer. In every confer- 
ence of this type, the usual wall plaque which reads ''THINK" (or 
even “THIMK”) should be replaced by a more important one (Figure 
1), “WHAT’S THE QUESTION?” The statistician's province is to 
tailor one of a stock of ready-made designs to fit the question, the 
laboratory worker makes his observations to carry out the design; 
but to set the question in terms which have laboratory answers is a 
problem which may take hours of work. Bicking*® has published a 
check list for test programs; the required laboratory work is step 
30 of 47 steps. 

Once the question has been agreed on, two quantitative values 
must be set before the statistician's work begins. The first of these 
is the reliability of the test method; the second is the least differ- 
ence in results which will be of practical value to the user. Also, 
this is the time to consider all the variables which might affect the 
results. Touchin's list? of the variables which affect exposure tests 
of paint panels is a good example; all these variables must be meas- 
ured or controlled. 


Preparation of the finishes. 

The panels. 

The human element in panel preparation. 
The conditions during painting. 

The conditions during drying. 

The aging period before exposure. 
Location of the exposure site. 
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8. Type of rack. 

9. Position of the panel on the rack. 

10. Weather immediately after exposure. 

11. Weather over the complete period of the test. 
12. Human element in evaluating breakdown. 





WHAT'S THE QUESTION? 


Figure 7 


Without good planning, it is quite possible that valuable differences 
among paints (in this case) will be completely obscured by weather 
changes, different observers, and the like. An expensive program 
will have been spoiled for lack of good planning. 

The outcome of a design conference should be the selection of a 
type of test program and the preparation of a detailed schedule of 
testing. The test schedule must be closely checked by the engineer, 
since he must attain in practice the experimental conditions which 
the statistician specified in theory. This does not mean that the 
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engineer should immediately reject as impractical all the unusual 
operating conditions which the test schedule demands, but a com- 
promise must be reached between the impossible and the impotent. 

Once the test schedule is set, it should not be changed in the 
slightest detail without the advice of the statistician. Small things 
are important—the interruption of work over a week-end or over- 
night; order differences in a routine operation—there will always be 
small points which were not covered in the plan. Most of these dif- 
ferences can be controlled and discounted if the statistician is aware 
of them beforehand, but some may be important enough to modify or 
destroy the whole program. 


EVALUATION OF RESULTS 


The biggest technical gap between the engineer and the statisti- 
cian lies in the interpretation of results and the presentation of data. 
Statistical terms are mathematical, exact, and usually unintel- 
ligible; their engineering interpretations, to be equally exact, always 
include the element of chance that is expressed by a probability 
calculation. But the interpretations must be made, and the engineer 
must make them. Moreover, this area of uncertainty covers the one 
part of engineering research in which statistics gives the most 
help—sure-thing probabilities usually need no interpretation. 

In a survey of statistical forms for the presentation of data, the 
most common one may be something of a surprise. It is a small 
thing, but its common adoption would double the amount of informa- 
tion available from laboratory work (Figure 8). It is this—that 
every average which represents a series of results should, as a 
matter of course, be accompanied by its standard deviation. To- 
gether, the average and the standard deviation are the best way to 
condense a series of measurements into two results. Alone, the 
average gives a false impression of accuracy; also, it is useless for 
evaluating the significance of differences. 

Significant differences are expressed as probability fractions. 
A probability fraction of, say, 0.95 indicates that chance variation 
alone would account for the difference found, in a sample of the size 
taken, only five times per hundred samples, if the sampling is re- 
peated indefinitely. Practically, this figure is used to imply that 
there is a 95 percent chance that the difference measured between 
two averages is due to a real effect of some changed condition, 
rather than to experimental error alone. 

Several parts of the x, y chart of a series of measurements have 
statistical interpretations (Figure 9). The correlation coefficient, 
r, is a number expressing the degree of fit of the points to the cal- 











DAVID N. COLLINS 








AVERAGE AND VARIANCE 


mere eee 
Figure 8 


culated line. When r = 1, the fit is perfect. If no line can be drawn 
through a cloud of compietely random points, then r=0. The r may 
be positive or negative without affecting its value, depending on 
whether the slope of the line is positive or negative. 

Both the slope and the intercept of the x, y line are statistical 
quantities. Each has its calculated standard deviation which can be 
used to make comparisons with other slopes and intercepts. For 
example, Youden* uses an x, y line to calculate results in chemical 
analysis. Here the slope (grams CaO per gram of sample) is the 
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CORRELATION 











X 


Figure 9 


desired result, and the intercept is the "blank" error, if any. A 
probability figure for the existence of a blank can be obtained by 
calculating the probability that the intercept is different from zero. 
The result of the analysis of variance is a table of variances, in 
which the effects of the rows, columns, and other subdivisions, and 
their interactions, are listed and compared with the experimental 
error. The table is packed full of information, but requires expla- 
nation to a user who is not familiar with this method of presentation. 


HOW TO DISBELIEVE A STATISTIC 


Statistics in research have reached the same degree of accept- 
ance in the 1950's as business statistics did about 1820. These 
produced a series of loud explosions from British politicians and 
businessmen, like Prime Minister Canning's. "Egad, I can prove 
anything by statistics except the truth!" But, like today's research- 
ers, the businessmen were forced to use statistics because those 
who used them were much better businessmen than those who did not. 
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Turning to research, the misuses and perversions of business 
statistics may also appear in the interpretation of research results, 
Darrell Huff, in a book called How to Lie with Statistics,® gives five 
qualifying questions for any data; these apply to corrosion results as 
well as to lung cancer. 

WHO SAYS SO? The laboratory with something to prove for the 
sake of a theory is dangerously liable to the suppression of unfavor- 
able data, or at least, to the presentation of the data in the most 
favorable light. Watch for the use of the ''O.K. Name," that is, the 
reference to the well-known author or well-known method which has 
been misapplied beyond its capacity. 

HOW DOES HE KNOW? Watch for evidence of a biased sample, 
one which does not fairly represent both sides of the question to be 
solved. This may not be conscious bias; for example, the choice of 
one panel from a pile never falls on the first or the last, but con- 
centrates near the middle. 

WHAT’S MISSING? Usually, what's missing is the number of trials 
on which the conclusion is based, or a reported difference in aver- 
ages without a measure of variance on which the difference can 
be judged. The use of percentages without the raw data is a good 
issue clouder, whether intentional or not. For example, there is 
the report that in one year 33percent of the women students at Johns 
Hopkins married faculty members. And so they had—one of the 
three students registered. 

DID HE CHANGE THE SUBJECT? This is the first of two questions 
which cover the engineer's problem in evaluating statistical results, 
Watch the validity of the step which moves from results to con- 
clusions. For instance, clocks A and B keep perfect time. When 
A shows nine o'clock, B strikes. Did A cause B to strike? These 
are two actions which are parallel, but which have no relation to 
each other. 

Finally, the second of the engineer's questions, DOES IT MAKE 
SENSE? Watch for extrapolations, either above or below the range 
of measurement. Watch for calculating errors. In expressing var- 
iance, the standard deviation associated with an average is calcu- 
lated to include two-thirds of the individual measurements within 
+1s.d. of the average. Smaller or larger values are usually in- 
correct, 


SALVAGE OPERATIONS 


Once the data is taken, statistical evaluation is of little help 
except to condense the raw figures into manageable form. In a large 
testing program there is usually enough variety in the conditions to 
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permit several statistical comparisons, and the benefit of statistical 
planning is in the amount of work saved. However, most tests of 
average size have a comparison missing, or do not reduce the ex- 
perimental error enough to prove real differences among the com- 
parisons. 

There are several methods of making the best of imperfect data. 
The most precise (and the least common) is to decide on the accept- 
able answer beforehand and select results to fit. Don't laugh; ex- 
amine your Own conscience. Without lecturing on professional 
ethics, I say this method is a complete waste of time and money. 

A second method might be expressed with the words ''Have 
Faith.'' Without enough data to prove the result, the conclusion is 
accepted as though it had been proved. This is the most common of 
all evaluations. 

The faith of the engineer in his accuracy is the biggest single 
deterrent to statistical methods. Methods based on faith save lab- 
oratory time and money; and if the results are easily foreseen there 
is no great danger of making an incorrect conclusion. Even if the 
results are somewhat confusing, the conclusions may be qualified by 
use of the verbs "indicate" or ''appear to indicate."' In any case, a 
little laboratory work is better than none and gives an air of re- 
spectability to an otherwise unfounded conclusion. 

There have been statistical salvage operations which produced 
good results. Both Bainbridge® and Lord” found their salvage so 
successful that the new conclusions were worth publishing for their 
technical value, as well as for a good example of statistical salvage. 

Not all salvage work will be concerned with mistaken or ill- 
defined conclusions. Unless they are closely checked, electrical 
calculators give wrong answers with just as much speed and finesse 
as they give the right ones. Mistakes are made in transcribing fig- 
ures from calculation to calculation. When conclusions depend only 
on numerical results, as these do, a small error can be deadly. 
Fyall’s study® of the determination of size on glass fabrics contains 
four calculating errors in his basic data, all of which should have 
been caught and corrected by cross checking before the analysis of 
variance was started. 


A LABORATORY UNDER CONTROL 


Engineers in general have their troubles, and the suggestion that 
they should be burdened with statistical techniques would not be 
made except for two real improvements which follow. First, sta- 
tistics save research time through elimination of the inconclusive 
project. Everyone can cite at least one such project—months or 
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years of data, progress reports which conclude, ''More data will be 
required before final conclusions can be made." With lessening 
interest, there is the final report which sweeps all the crumbs of 
information into a pile, saying, ''The data appear to indicate... ," 
The data indicate nothing much, really. Ten thousand dollars spent 
for no conclusions. 

The second benefit from statistical methods is that variance 
measurements of standard laboratory techniques can be calculated 
and reported and can be used for planning later work. This means 
smaller test plans, or increased precision, or both, as each test 
draws on previous results to prove significant differences. In a 
laboratory which uses all its data, the variance is as important as 
the average with which it is associated, and the report of an un- 
usually high or low variance should raise as many eyebrows as the 
report that three Super Gas Savers in series will save 125 percent 
of your gasoline. 


WHERE TO GET HELP 


Sucha revolution in laboratory operation will not come overnight, 
and will not come at all without expert advice. The best expert ad- 
vice for day-to-day operations can be applied through the use ofa 
statistical advisor at the lowest level of research planning. This 
advisor should not be a mathematician, but should be an engineer 
with statistical training who can apply his full time to advising on 
plans and to interpreting data. This is the’'plant engineer" concept 
applied to research work—the man responsible for patching leaks 
and solving problems to keep up production. One man like this, plus 
a calculator operator, can handle the problems of 50-75 research 
workers. 

The more-senior research planners who cannot spare time for 
training in the techniques may have to accept the techniques on faith 
for the results they bring. The technical literature of the last ten 
years is sprinkled with expository articles on ''the new statistical 
approach,"’ but none of them (including this one) have space enough 
to do more than suggest that statistics is a good thing and to illus- 
trate some of the techniques. Do-it-yourself from books is slow, 
there is still no royal road to learning. However, most of the books 
on statistical methods will repay even a few minutes taken to read 
the Preface and Chapter I. Beyond this, the best results are ob- 
tained by classroom training. 

Some training in statistical methods should be part of all scien- 
tific education. The half-year mathematics course which is the in- 
troduction to statistics seems quite satisfactory for teaching basic 
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techniques, but any application to the problems of laboratory work 
would be out of the question. This application, for chemists, should 
come quite casually in the undergraduate courses of quantitative 
analysis and physical chemistry, both in laboratory and classroom 
work. Some attention to the problems of variance in laboratory ex- 
periments would also help make laboratory training less formalistic 
and more practical. 

For the graduate training of those who plan to specialize in the 
engineering applications of statistics, mathematics courses are 
already in use, built around the statistics of business. However, 
there should also be a graduate course in research planning to fit 
the special nature of laboratory statistics. 

As more engineers are formally trained in statistics, there will 
be less of a gap between the engineer and the statistician. However, 
the demands on statistical specialists will increase, as present 
plans become more complex and new ones are developed to better 
fit the requirements of engineering research. As trained men be- 
come available, those of us can be retired who (like myself) have 
little talent for mathematics, but who are trying to fill the gap until 
better men arrive. 
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Determination of the Necessary Sample Size in 
Life-Testing When the Length of Life is 
Exponentially Distributed 


By Calvin J. Kirchen, Engineering Staff 
Ford Motor Company, Dearborn 


FOREWORD 


The necessary information for the determination of sample size 
for experimental or inspection work is given. Its use is exempli- 
fied in sampling inspection as an orientation. The major part of the 
paper is devoted to its application in life-testing work when times 
to failure have an exponential distribution. Formulas for the cal- 
culation of various quantities needed in this work and tables of im- 
portant quantities are given. 


At the beginning, it is necessary to state the basic concepts in- 
volved in sample-size problems. They will be stated in terms of 
the sampling inspection operation that is almost universally used 
in the receiving departments of our many industrial plants. 

Suppose a shipment of 20,000 bolts is to be inspected. Since 
only a sample is to be inspected, it is not only possible but quite 
probable that there will be defective bolts which will not be nicked 
up in the sample. In the event the shipment is accepted by sampling 
inspection, these defectives will remain in the shipment, to be dis- 
carded perhaps on the assembly line. The point of this is: requir- 
ing a shipment of 20, 000 bolts to be 0 per cent defective is actually 
unrealistic, because not even 100 per cent inspection could guarantee 
it. Hence, some small per cent defective should be defined as “ ac- 
ceptable” per cent defective; this per cent defective is usually re- 
ferred to as a “quality level,” in particular, as the “acceptable 
quality level”: AQL. Obviously, there must also be a specification 
of an “unacceptable quality level”: UQL, and, of course, the UQL 
in per cent is greater than the AQL in per cent. It does not follow 
incidentally, that UQL is the minimum feasible per cent defective 
greater than AQL. 

A little reflection on the fact that the sample may contain more 
than the expected number of defectives and hence indicate that the 
shipment of bolts being sampled is of unacceptable quality (an ex- 
ample of “sampling error”), when it is in fact acceptable, will make 
reasonable the concept that in drawing conclusions from sample 
data there is the unavoidable risk of rejecting an actually acceptable 
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shipment. This risk is always denoted by the Greek letter alpha, 
a, and is specified at the discretion of those administering the in- 
spection operation It is customarily given as a percentage and 
rarely exceeds 10 per cent. It is important to associate AQL anda 
closely; in fact, they are inseparable. Similarly, there is associ- 
ated with the UQL the risk of accepting a shipment of unacceptable 
quality; it is denoted by beta, 8, and is likewise specified by the 
inspecting organization, with a value usually not exceeding 10 per 
cent. It may or may not be equaltoa. A widely adopted system 
fixes a at 5 per cent and £8 at 10 per cent, but their values are arbi- 
trary. In statistical language a is usually termed a Type 1 error 
and 6 a Type 2error. In inspection circles they are known respec- 
tively as producer’s risk and consumer’s risk. 

The relations of these four concepts are clearly shown by a graph 
which gives the probability of acceptance by a sampling plan of ship- 
ments of various grades of per cent defective. Such a curve is 
called the operating characteristic of the inspection plan. The op- 
erating characteristic of any inspection or experimental procedure 
for specified values of AQL, a, UQL and 8 can and should be calcu- 
lated and graphed. 

If we set up the following specifications for a sampling inspection 
plan: AQL = 2 per cent, UQL = 9 per cent, a = .05 and 8 = . 10, then 
standard statistical methods or a reference to F. L. Grubbs, “On 
Designing Single Sampling Inspection Plans,” Annals of the Institute 
of Mathematical Statistics, 20, No. 2, 242-256, will show that in- 
specting a sample of 70 pieces taken at random from a lot and ac- 
cepting if no more than 3 defective pieces are found, will satisfy the 
four requirements above. The operating characteristic is shown in 
Figure 1. 

The remainder of this paper will be devoted to the use of these 
concepts in life testing, with particular emphasis on the case of an 
exponential distribution for the times to failure of items tested: 


le-* L, 


L 


where x is measurement of life on some scale (units of time, or 
number of test-cycles) and L is average life of the material or kind 
of item tested. The methods described are based upon the paper 
“Life Testing’ by Professors Epstein and Sobel, Wayne State Uni- 
versity, in the Journal of the American Statistical Association, 48, 
No. 263 (1953), 486-502. 

The basic concepts in terms of life testing are: 


1. The minimum life (number of cycles of testing or life on time 
scale) a part must last to be acceptable: La. Associated 
risk of rejection isa. 
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Figure 1 


Operating Characteristic 


2. The maximum life a part may last in being unacceptable: 
Ly: Associated risk of acceptance is 8. 


The notation given above (La and Ly) is different from that used 
by Epstein and Sobel, who used notation customary in mathematical 
statistics, which is frequently less suggestive of actual applications. 

While it must be emphasized that an exponential distribution of 
times to failure is basic to the methods to be discussed, neither 
space nor time is available to discuss this further. See D. J. Davis, 
“An Analysis of Some Failure Data,” Journal of the American Sta- 
tistical Association, 47, No. 258 (1952), 113-150, with particular- 
ly the middle graph of the bottom three on p. 116 and the last para- 
graph on p. 117. 

Professors Epstein and Sobel had particularly in mind the type 
of testing common with electrontubes where all items of the sample 
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are tested simultaneously so that the first tube to fail marks itself 
unmistakably as the weakest of the sample; the second tube to fail 
as the second weakest, and so forth. Under such conditions it is 
not necessary to test until all items have failed, to reach a decision 
as to acceptability or unacceptability. Using the notation xj_» for 
the life of the i-th item in a sample of n items, Epstein and Sobel 
showed that the statistic calculated from the life-data of r failed 
items out of n under test, which maximizes the exponential distribu- 
tion given above, is: 


Xi,n + Xan +eecvs + Xr in + (n-r)Xr.n 
Lin = “ 





It is apparent that this is similar to an average lifetime of the r 
failed items alone, but it differs from that because it is weighted 
by the final term. It can be remarked here that the calculation of 
L », 48 given above is the only calculation that need be made with the 
Xi,n in order to reach a decision. An unusual feature of the testing 
procedure that should be emphasized here is that r, the number of 
failures, is chosen in advance of testing and n, the number of items 
under test, is left arbitrary. 

The natural question, “What guide, if any, is there for the choice 
of n, the number of items to be put under simultaneous testing?’ is 
answered by Table 1 of Epstein and Sobel’s paper. For various 
values of r and n, this table gives the average amount of time it 
takes to attain the rth failure in n compared with the rth failure in 
r. This table is calculated from: 





i=l 


An excerpt from Table 1 appears as follows: 





r\n_ 1 2 3 4 5 10 15 
Bev: bed . 50 . 33 . 25 0 .e et 
Bg ee ee . 56 . 39 28 Mb: 1 
$e - 1. 00 . 59 ZS «<*> oe 
pe - - 1. 00 oe 
5 1.00 .28 
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Assuming a rational basis for the choice of r, the experimenter 
has here a guide to making a choice from the many available values 
ofn. Time-wise, it will be advantageous to make n somewhat larg- 
er than r; this will have to be balanced against the cost of obtaining 
items for test and for testing them. 

The final phase of this work is the definition of the test of Ly, 
whereby a decision can be made as to the acceptability of the items 
under test. Explicitly stated, what is desired is how “best” to use 
the first r ordered observations of a sample of n items in order to 
discriminate between L, and Ly, La >Ly. “Best” has the meaning 
in statistical testing of hypotheses of specifying a test which has the 
property, among all tests having a fixed probability a of rejecting 
L=L, when true, of having the largest probability of rejecting 
L = La when L = L pis true. 

Methods of mathematical statistics are employed to obtain a test 
value denoted by C, which is expressed in terms of a well-known 
statistical distribution called chi-square, which is tabulated in most 
statistical texts, the acceptable life La, and the quantity r: 


C = Lag Xi-a (2r)/2r. 


For acceptance of items under test, we must have L,, > C. 

The factor which resembles “x square” is “chi-square.” Both 
the subscript (l-a) and the 2r in parentheses require comment here. 
The shift from a to l-@ was made for convenience in the mathe- 
matical statistical derivations. What is more important to the user 
is the fact that there is at least one set of tables of the chi-square 
distribution which defines the probability levels as complements of 
those Epstein and Sobel had in mind. This set is in A. Hald, Sta- 
tistical Tables and Formulas, New York: John Wiley and Sons, Inc., 
which supplements their excellent work Statistical Theory with Engi- 
neering Applications (same publishers). After remarks on the 
meaning of (2r), a means of identifying the proper set of tables to 
use will be given. Chi-square tables are arranged in columns, one 
of which is headed by n (extreme left column) and the remainder by 
probability values. The entries in the n column are known as “de- 
grees of freedom,” and the quantity (2r) defines the number of de- 
grees of freedom for which chi-square is evaluated. The number 
of components to be tested to destruction is, then, half the degrees 
of freedom. Since .05 is a frequent probability involved, it will be 
used to distinguish between Hald’s table and those Epstein and Sobel 
had in mind; if the .05 probability column for 1 degree of freedom 
has a chi-square value of 3.84, then a = .05 in Hald's sense. 
In Hald’s chi-square tables, 3.84 appears in the .95 probability 
column for 1 degree of freedom. 
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For reasons of practicality both in theory and application, the 
ratio L,/Lyor its reciprocal, and the value of La, are used rather 
than specific values of La and Ly. Values of r which are used are 
smallest integers such that 


Xia (2r)/*3 (2r) > Ly/La e 
From the definition of the critical value of C above, we have 
C/L, = X?_q (2r)/2r. 


Hence C can be evaluated in advance of experimentation. If L,,, 
calculated from the experimental data, exceeds C, then we con- 
clude the material is acceptable, running the risk 8 of accepting 
material whose life is Ly. Epstein and Sobel made up nine sub- 
tables, using all combinations of .01, .05 and .10 for both a and 8, 
and with La/Ly specified for the following series: 1.5, 2, 2.5, 3, 
4, 5, and 10. In each sub-table can be read the appropriate value 
of r and C/L,. Their table fora = .05, 8 = . 05 follows: 


L/L, iF C/L, 
15 67 8079 
2 23. 6834 
2.5 14. 6046 
3 10 =. 5426 
4 7  .4694 
5 5 .3940 
10 S <a 


As an example of its use, if we assume La = 100 and Ly = 20, then 
La/Ly= 5 and hence r= 5. C = La(. 3940) = 100 (. 3940) = 39. 40. 
If we choose n = 10, then from Table 1 we see that 5 failures in 10 
tested would be obtained on the average in . 28 of the time required 
to observe 5 failures in 5 tested. This time ratio can be reduced 
even more if more than 10 are tested. At any rate, acceptance 
depends on L,, > 39.40. To illustrate the method of obtaining r 
and C for other than the given values of La/Ly, assume we want 
this ratio to be 8, witha =.05, 6=.05, soLy La=.125< Xi-g 
(2r)/X% (2r). Comparison of the l-a@ = 1-.05 = . 95 column and the 


. 05 column of achi-square table shows that n = 6 is the least integer 
such that .125 < 1.64/12.6, so 2r=n = 6, or r = 3, and C/La = 
1.64/6 = . 27. 

As a final remark, the calculation of the operating charactertistic 
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of a life-testing experiment is shown by an example. The procedure 
is indicated by relation (27) in Epstein and Sobel's paper. 
Probability of accepting L = La when L is true value = 
La 


Probability Ix? (2r) > i: Xi-a (on) ‘ 


Following the example of the operating characteristic shown for 
the sampling inspection plan, Figure 1, the probability of acceptance 
is plotted as the ordinate, with the ratio La/L as the abscissa. Let 
a = .05, and a specific value of r be ehosen; then Xj-q (2r) is a 
specific value in the chi-square table. Choose a value of L and ob- 
tain La/L. Multiply this by Xj-q (2r) = X75 (2r). Find the value of 


x*(2r) equal to x". 9s(2r). Note the probability level this value 


of X* is in. Plot the probability level for the value of La/L. It 
will be found quickest to obtain a series of values for La/L by cal- 
culating x*(2ry/x*, os (2r) for the various values of X*(2r) > X” os (2r). 
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The Use of Computing Machines 


By Claus B. Ludwig 

Motor Truck Engineering Department 
International Harvester Company 
Fort Wayne, Indiana 


INTRODUCTION 


One and a half centuries ago, the substitution for manual effort 
of heat energy for power plants inaugurated the steam age. Then 
nobody could foresee the incredible changes in business, production 
and machinery that would come into being as a result of the so-called 
industrial revolution. Today, industry is again experiencing funda- 
mental changes which can be considered a new revolution. In order 
to meet present-day requirements of cost reduction and increased 
productivity, automatic control and data processing are now being 
employed in manufacturing and clerical operations. This trend 
toward increasing efficiency is growing with the improvement of 
machines which can perform the functions of control and data pro- 
cessing. An industry which does not recognize this trend will cease 
to exist. 

Although self-control devices have been used for some time, only 
the discovery and fast development of electronics and its general 
application to feed-back systems have made it possible to improve 
control equipment enough to stimulate revolutionary changes in in- 
dustry. The great evolution in the field of data processing is also 
due to electronics. During the war, new kinds of electronic com- 
puters were developed wherein the mechanical computing parts were 
replaced by magnetic and electrostatic units. The significance of 
these computers lay in their great flexibility and exceedingly high 
speed. It became possible to solve extensive mathematical prob- 
lems in an unbelievably short time and to process tremendous 
amounts of statistical data. We know today that these electronic 
computers were used to solve such complex problems as those per- 
taining to the path of guided missiles and of projectiles for aviation 
defense, weather forecasting based on data from many stations, and 
atomic power development. The presence of complex problems in 
industry has also animated the application of the electronic com- 
puters in peacetime work. 

This paper pertains to the applicationof automatic computers for 
design analysis and for the prevention of failures of machine parts. 
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It is, therefore, necessary to go into detail about automatic com- 
puters. In doing so, I have to talk about one of the most fascinating 
inventions made by man. 

We all know of the two extremes public thought takes concerning 
these new electronic computers. One is that of disregarding most 
of the merits of the computers and reproaching the inventors for 
claiming their machines can “think.” The other is that of waxing 
enthusiasm which claims that no problems will remain unsolved in 
the future and that “giant electronic brains” can do anything. 

As always in life, the truth is somewhere in the middle. These 
machines have a great many merits, and, as we will see, their way 
of doing things can very well be regarded as thinking. On the other 
hand, these machines are far from being a “brain,” although, the 
function of a brain is also thinking. To clarify, we can categorize 
the process of thinking into original thinking and routine thinking. 
Original thinking might make up, say, only 10 per cent of our mental 
effort and routine thinking the remaining 90. The general purpose 
of the electronic computers is to relieve men of the burden of the 
90 per cent routine mental effort, but men must themselves still 
supply the important 10 per cent original mental effort to put the 
computers into use. 

We engineers are now faced with a great variety of small and 
large-scale electronic computers and with the yet greater task of 
appraising them. We would like to know what these computers can 
do for us, what their limitations are, what effort we have to make 
in order to apply them most efficiently, and so forth. I want to 
discuss these points and to present an outline of engineering prob- 
lems which may be solved by electronic computers. 


PRESENT SITUATION IN ENGINEERING 


The load and speed requirements of machine parts have _ in- 
creased tremendously in the past fifty years, while the size of 
machine parts has decreased. This development could not have re- 
sulted from sheer guesswork. It came about not only through prac- 
tical experience but also through analytical work when it was 
impossible to find the optimum design by trial and error methods. 
Analytical work may be defined as the statement in mathematical 
terms of the behavior of material and related systems. The quality 
of this analytical approach depends upon the assumptions which are 
introduced into the theory. On the one hand, it is desirable to ap- 
proach the exact theory as closely as possible by using accurate 
assumptions and exact boundary conditions. On the other hand, it 
is necessary to keep the assumptions and boundary conditions as 
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simple as possible in order to allow for a numerical solution of the 
problem. Therefore, the analytical approach is a compromise be- 
tween the best theory available and the simplest theory possible. 

In the first approach to a given problem, the simplified theory 
will outweigh the more exact theory and the answer may be equiva- 
lent only to a rough estimate. Then, as time and people improve 
the first approach, the answer becomes more and more refined until 
sooner or later the limit is passed beyond which the theory does 
not provide an economical answer, i.e., an answer which can be 
evaluated numerically in a reasonable time. 

In the analytical approach, the language of mathematics is used 
to condense the entire concept of the investigation into a few short 
symbolic statements, but in the thorough investigation lengthy cal- 
culations may involve an enormous quantity of numerical data. I 
think it is safe to say that the mathematical solution of a given prob- 
lem may require 10per cent of mental effort in the sense of original 
thought and that the numerical evaluation of this solution may re- 
quire the balance of 90 per cent of routine thought. Men will always 
be happy to spend the 10 per cent of original thought but the burden 
of 90 per cent routine thought kills nearly all of the more complex 
projects. 

Now, in this situation, new helping tools are given to us which 
eliminate the monotonous and time-consuming computations and 
which can be the essential bridge between an engineering concept 
and its numerical solution. These new tools are the automatic 
computers which are completely controlled during operations of 
data processing by the use of stored programs in electrostatic 
memory or by the use of wired control panels. 


SURVEY OF APPLICATION OF AUTOMATIC COMPUTERS 


Electronic computers may be classified as either digital or 
analog computers. A digital computer is a device which operates 
with discrete values given in the form of digits. The accuracy of 
the performed result depends upon the number of digits the machine 
can handle. With the large scale digital computer, the precision 
can go up to approximately 150 digits. An analog computer can 
handle only physically measurable quantities. Therefore, the dis- 
crete values which are used in a digital computer have to be con- 
verted into quantities such as lengths, currents or angles for the 
purpose of calculating with an analog computer. 

The difference between a digital and an analog computer can be 
explained by comparing a desk calculator and a slide rule, respec- 
tively. The slide rule is an analog computer because the logarithm 
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of numbers is represented by lengths. A desk calculator is a digital 
device because it can handle only discrete values in the form of 
digits. For both the analog and digital computer, the problem must 
be stated before it can be programmed. The difference is that the 
answer does not have to be known for an analog and the approximate 
answer must be known for a digital computer. 


Analog Computer 


In general, the analog computer can be applied in basic research 
studies. It will be used more for the investigation of the behavior 
of physical systems than for the evaluation of specific design data, 
since it has a limited accuracy. However, in problems where an 
accuracy with four digits is sufficient, an analog computer might 
be applied very efficiently. Of course, time will be required to 
translate the given information from the form of digital data into 
that of physical quantities. Also, the quality (i.e. the accuracy) of 
this translation depends upon the skill of the programmer. Figure 1 
gives an outline of the applications of engineering problems oc- 
curring especially in the automotive industry soluble with an analog 
computer. 


Figure 1 


Analog Applications* 


1. Parameter studies, that is, analysis based on insufficient data or where we 
must design by trial and error; the method involves the systematic sub- 
stitution of different combinations of parameters which will be the most 
useful for achieving a desired result. 

2. Vibration of structures. (For example torsional vibration in crankshaft.) 

3. Spring systems and vibration absorbers. 

4. Dynamics of mechanisms 

5. Fluid flow problems. 

6. Analytical studies involving gears or other mechanical linkages requires 

the determination of backlash. Backlash angles from zero to any practical 

value are easily introduced to the system. 

Engine crankshaft main bearing loads. 

8. Lubrication of engines, where the problem is to determine the quantity and 
pressure of the lubricating oil at various points throughout the lubricating 
system 

9. Solution of sets of differential equations. 

10. Gas flow in a fuel system and fuel combustion. 

ll. Behavior of valve in an engine. 

12. Performance characteristics of engineering systems. 

13. Simulation of thermodynamic cycles 

14. In the farm implement field, the design of machine parts, such as knife 
blades, to give the desired trajectory of the material being handled. 


ompiled by Miss H. 8S. Siegel of the International Harvester Company 
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Digital Computer 


In general, the digital computer can be applied for the numerical 
evaluation of mathematical formulae. Since the digital computer can 
execute only the four basic arithmetic operations of addition, sub- 
traction, multiplication and division, all other operations, such as 
integrations, root extractions, differentations, evaluations of angles 
from trigonometric functions, have to be expressed in terms of these 
four basic operations. Therefore, considerable time is required to 
program given formulae. But we will see that, with a digital com- 
puter, the evaluationof a sine by using fifty terms of a powerseries, 
for instance, requires actually much less time than to look up this 
value in a table. Furthermore, the automatic computer eliminates 
the errors possible with the use of tables. 

Generally speaking, there are two types of computations involved 
in technical problems: 


1. Putting numbers into formulae which are either exact or close 
approximations. 

2. Trial and error method when explicit formulae are replaced 
by a set of formulae and the calculations have to be repeated until 
the result is found. 


Figure 2 gives an outline of the problems soluble with a digital 
computer by EDPM because of their monotonous part of numerical 
calculation. They are arranged into two groups as follows: 


1. One group includes current problems which are handled at the 
present time by a desk calculaior. 

2. The other group consists of problems which were never at- 
tempted because of the lack of ability or time. 


Obviously, this outline is meant only to be suggestive. Only ex- 
perience will give a more complete list of all the possible engineer- 
ing problems. 


ELECTRONIC DATA PROCESSING MACHINES 


From here on my paper will deal only with large-scale digital 
computers (EDPM) which can also be used as data processing ma- 
chines. Examples of these automatic computers are the Univac by 
Remington Rand, CPC, 702 or 705 by IBM, Mark IV by Harvard 
University, Enivac and Edvac by Moore School of Engineering, and 
some others. These machines can be purchased or rented at dif- 
ferent costs which may vary on the rental basis from $3,000 to 
$70,000 a month. Another possibility is to use a large-scale com- 
puter in one of the central computing laboratories at a rate of ap- 
proximately $500 per hour. 
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Figure 2 
Problems soluble with a digital computer by EDPM 
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The question arises: What type at which price shall anyone 
choose? This, of course, depends upon the size of the engineering 
outfit and the complexity of the problems which occur and also on 
whether or not the computer will be used for accounting purposes as 
material requirements andinventory control, maintenance problems, 
payroll jobs, etc. Careful consideration has to be given to all 
these factors before renting or purchasing one of these machines. 
General practice has been to send people to the schools conducted 
by the manufacturers of the computing machines to ascertain the 
feasibility of the machines for specific problems. 


Short Description 


As an example of a computing machine, the IBM 702 is shown in 
Figure 3. 

Imput Units. Data are given into the machine from both magnetic 
tape and punch cards. 

Output Units. Results can be written on magnetic tape, printed 
on paper, punched into cards, or plotted as curves from a television 
tube. 
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Figure 3 


An Installation of IBM Electronix Data-Processing Machines 
Type 702 and Associated Equipment 
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Central Units. All mathematical operations are performed in H 
the central unit which contains the electrostatic memory and the | solve 
accumulators and the arithmetic and logical unit. exam 

Drum Storage. Additional storage is possible with the use of becal 
magnetic drums. will « 

Operations. Besides the four basic mathematical operations purp' 
(addition, subtraction, multiplication, division), many logical op- swer 
erations can be executed by the machine. These include the modi- of th 
fication, under given circumstances, of the instructions and de- | T 
cisions on such questions as the following: Are two figures equal, is th 
higher or lower? Is one figure minus, plus or zero ’?and Is the re- the h 
sult right? macl 

Internal Checking Procedures. Many possibilities of checking lem 
the internal operations are provided which give a control of ac- lems 
curacy of results. , «Seri 

Stored Program. The instructions included in a program con- the 7 
tain five positions. The first position is the code of the operation woul 
which has to be performed. The remaining four positions represent of th 
the memory address at which the mathematical data is stored. shall 

Since the program is stored in memory in the same way as data, give 
the machine may modify instructions by substituting a new address. mac! 
Therefore, a program can operate on itself and compute one ofits | I 
own instructions. ther 

Time of Operations. As an example of the high speed of the prob 
large-scale digital computers, the time required for some of the prob 
operations is of interest as it is given in the manual for IBM 702: prob 

shou 


resu 

237, 000 additions or subtractions 

(5 digits) per minute 

' Prep 

50, 000 multiplications (5 digits x 
5 digits) per minute ' I 
' prob 
24, 000 divisions (5 digits) per minute | quir 
434,000 logical decisions per minute a 
' orde 
The Limitations of the Large-scale Computer be F 
| Figu 
The limitation of the digital computers is the human brain. The | gine 
machine willnever be more clever than the human brain which works digit 
with this machine. The machine will never be able to give a result grar 
which could not be reached by the human brain also. That is the ; (i. e. 
limitation: What the human brain cannot figure out, the machine mate 


cannot either. 


SP 
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However, in many cases, the human brain would not be able to 
solve a problem in the time in which the answer is required. For 
example,the machine will not solve one equation with two unknowns 
because the human brain could not do it either. But the machine 
will calculate in one second the path of an airplane for air defense 
purposes where an individual would need several hours. Both an- 
swers will be right, but the human answer will be useless because 
of the time factor. 

There is another interpretation of the statement, the limitation 
is the human brain. Not only is the machine no more clever than 
the human brain, but the human brain must be more clever than the 
machine. The human brain must have adefinite concept of the prob- 
lem which the machine shall solve. For example, one of the prob- 
lems in our cam calculation is to find the most favorable power 
series. It is easy to make a program for the machine to calculate 
the 780,000 series. The machine would run about 1,000 hours and 
would give approximately 3 1/2 billion values. Without doubt, some 
of these values will be very favorable for our purposes, but who 
shall look for them? The great danger is that the machine could 
give too many values if there is no concept. These 1,000 hours 
machine time would be useless. 

In order to make an intelligent use of the “electronic brain,” 
there should be a definite concept of calculation for every single 
problem. This involves detailed theoretical investigations of the 
problems as well as programming. A program for one of these 
problems should never be based entirely on existing techniques but 
should be improved before programming for the EDPM. This will 
result in a better design and eliminate expensive test runs. 


Preparation of Programs 


In order to program a given problem for a digital computer, the 
problem must be stated in the language of the machine. This re- 
quires an exact knowledge of the whole problem. Even the approxi- 
mate answer must be known. All equations and formulae have to 
be broken down into the four basic arithmetic operations. Then, in 
order to facilitate the programming, a detailed blockdiagram should 
be plotted, including all the steps of the machine operation. See 
Figure 4. The work up to this point has to be done by a skilled en- 
gineer who is well acquainted with the problem itself and with the 
digital computer as well. Anybody sufficiently familiar with pro- 
gramming can transfer the block diagram into actual program steps 
(i.e., the machine language). An average problem with approxi- 
mately 1,000 instructions requires about two to four weeks for stat- 
ing the problem and plotting the block diagram and about four days 
for writing the actual program. 
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Next, the written program is transferred onto punch cards (or 
whatever is needed for feeding into the machine). These punch cards 
can be stored to be used later when problems of the same nature 
arise, provided that the factors relating only to a specific design 
are put on separate cards (see Figure 5). Therefore, in writing the 
program, it is wise to plan for all future possibilities, i.e., to 
make the proper spacing. Nevertheless, the great flexibility even 
allows for future alteration in the basic program. 

Another task is to prepare subroutines which include procedures 
and methods to be used over again in different problems. Examples 
of subroutines are extracting roots, obtaining the angles when the 
sines, cosines, etc., are given and vice versa, solving polynomials 
of any degree, etc. 


The Staff as an Essential Supplement to the Machine 


In order to fully utilize the large-scale digital computer for en- 
gineering projects, a definite program should be started in which a 
group of mathematician engineers become familiar enough with the 
potential of the equipment so that engineering problems may be ap- 
plied effectively. This group has to accomplish two purposes: 


1. To investigate all standard design procedures in which 
lengthy calculations are involved and to set up a new improved 
procedure. Many standard design problems can then be completely 
programmed and stored. If a new design is required, only one or 
two punch cards containing the specific design data and the code 
number of this program have to be transmitted to the machine. 
After a short time, the printed results will come back. The pro- 
gram remains in the storage room so that it is available for the 
next application. In this way, many design problems can be estab- 
lished within one hour while otherwise weeks of desk calculation 
would be required. 

2. To investigate problems which cannot be solved economically 
at the present time. For these problems, in most cases much time 
is needed to find the most favorable approach. 


Finally, it should be pointed out that the quality of the computer 
group determines the efficiency of the machine. 

There appears to be a definite advantage in having a central 
group of engineers with the greatest possible competence instead 
of having several engineers from various sections for the following 
reasons: 


1. The central handling of current design problems involving 
lengthy calculations will avoid overlapping. 
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2. Only a central group can coordinate the problems in such a 
manner that a great number of parts of EDPM programs can be re- 
used for new problems. 

3. The personnel in this group are kept in training for EDPM 
programming. 

4. Because of their education, they will be able to increase the 
theoretical understanding of most of the given problems. 

5. The people of this group will gain enough experience after a 
certain time so that they can apply the EDPM to unsolved problems 
most efficiently. 

6. Because of their permanent training in analytical investiga- 
tions, they can easily exploit their experiences in fields in which 
they were not previously familiar. They have the advantage of 
bringing a fresh approach to old problems; and, therefore, they 
might be able to find the general solution where only a small part 
was known before. 


EXAMPLES OF APPLICATION 


In this section, I will give a short survey of problems which have 
already been solved in the Motor Truck Engineering Department of 
the International Harvester Company. 

First, it might be interesting to compare the actual machine 
time required by the different computers for one specific problem. 
Take, for example, the cam profile calculation based on the poly- 
dyne method. In this method, the evaluation of 


Yo = hy +9, (h+%C,0) + $2 n(n-1)C,0"* 


has to be made where hr,h, 9, @ and Cn are the design data. gis the 
angle of duration and n indicates the powers p,q,r,s,t. There are 
approximately (all operations are executed with 7 digits) 


3,000 additions and subtractions 
4,000 multiplications 

800 divisions and 
1,000 logical decisions 


involved. The machine time to execute these operations was on the 


Desk Calculator 100 hours 


IBM 604 12 hours 
CPC 4 hours 
IBM 702 20 seconds. 


It appears that the time required for execution of this program 
for the different machines decreased tremendously as we used 
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more competent machines. This improvement is not only due to the 
higher internal speed but also to the ability of the larger machines 
to make logical decisions. 

I do not have available the exact figures of rental cost, but I do 
know that the IBM 702 is much more expensive than the CPC. There- 
fore, a question arises as to whether or not the time saving of 
roughly four hours is worth the additional money for renting the 
IBM 702. The decision is in favor of the larger IBM 702 for the 
following reason. There are many programs which require more 
time than the twenty seconds on the IBM 702. If we assume a linear 
relationship between the machine times for the two computers, al- 
though actually the CPC will slowdown for more complex programs, 
we have atime table. To illustrate, I will choose the evaluation of 
the entire solution of the differential equation which describes the 
valve gear train at non-design speed: 


d’y 2 2 n BY n-2 
— + py = BiA+ FCpO"| +(—) 2 n(n-1)Cn0 
d@’ n Bo n 


This solution was obtained on the IBM 702 after about fifteen min- 
utes. Assuming a linear relationship, the CPC would have needed 
180 hours. It appears that for more complex problems, the smaller 
computers are no longer desirable for engineering applications. 

Let me describe briefly the answer obtained by the IBM 702 in 
the case of the solution of the differential equation (the formula is 
shown in Figure 6). There were executed approximately: 


60,000 additions and subtractions 
120,000 multiplications 

20,000 divisions 

50, 000 . logical decisions 

10,000 figures printed with 7 digits each. 


The printed problem solution is shown in Figure 7. The answer 
consists ofthe evaluation of eighteen basic series of p,q,r,s,t 
which are needed for the correct selection of the power series and 
of the evaluation of the amplitudes cf the transient vibrations in the 
valve gear system for a great variety of power series and fora 
wide range of speed ratios. This work was necessary to determine 
the relationship between the selected cam profiles and the transient 
vibrations in the system in order to avoid valve failures. 

In the meantime, many more problems have been analyzed and 
programmed. Some of them required only the translation of exist- 
ing equations into machine language, others required an extensive 
analysis for several months. Many mathematical procedures and 
so-called “subroutines” have been employed for the problem 
solutions. 


ed 
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Step 1 
y= hy + by +Py” 


on , ‘ 
y = h(1+ 0N,x ) + VIN}x" + ADNnx" 
y = hSaNpx™ + VlapN px" + ALapNax” 


n(n-1) 
x =0/c; a = x 10° 





no = 7PaFSt/n ye v_ -E(n-6-w)100 





n 
m(m-n)/ a») m(m-n)/ wn) 


a 1000 
ny? ,——__—*— 
a m(m-n)/ wen) 


oO 


Vv a 
Jv/jpax = 142N2 +72 Nn + px Nn 
n+i n+vi 








Step 2 


S. 2m/92™ ( -@/2+.... (-1)™ &*”- Y2m-2)) 


S, = 2m-1/9" (6 ic eal (-1)™"? &™- ¥/(2m-3)) 


Figure 6 


Equations of Cam Calculation 


In Figure 8, a list of engineering automotive problems is given 
all of which have been successfully solved by the IBM 705. As far 
as possible, approximate time figures are also given. 

In the following, some basic ideas about subroutines are pre- 
sented. The significance of subroutines is based on the fact that 
they might be permanently stored on cards or tape and only the 
transfer addresses be known. They can then be used in the solution 
of any program involving subroutines. Subroutines might be pref- 
erably programmed in floating decimal point. But it is not always 
favorable to use floating decimal point mathematics in engineering 
problems. Therefore, either a conversion routine or a second sub- 
routine procedure with fixed decimal point should be available. 
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The Printed Problem Solution of the Differential Equation in Figure 6 
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Desk Calc. IBM 705 
Project Time Time 
Hours Minutes 
General Design Studies 
and Specific Design Solutions 
Torque Converter Design 40 1 
Performance Studies 8-40 2-1 
Cam Design 100 1.5 
Pressure - Volume Diagram 25 5 
Rear Propeller Shaft Data 3 .06 
Heat Regenerator for Gas Turbines ? 
Combustion Chamber Design 120 2 
Fuel Economy 
Gear Designs (Spur, Helical, Bevel, Hypoids) 
Vibration Problems 
Truck Vibrations 40 l 
Tractor-Trailer Vibrations (3-mass system) 50 l 
Transient Vibrations in Valve Train 25 5 
Engine Balancing ? 
Stress Problems 

Bearing Loads 250 6 
Spring Design 8 . 
Strength of Front & Rear Axle Shafts 8 
Bending Stresses in Frames 16 5 
Bending Stresses in Discs 40 1 
Stresses in Rotating Discs of 

Non-uniform Thickness 40 l 
Stresses in Gear Teeth (bending and 

compression) 100 2 

Statistics 
Statistical Analysis of Data 4-40 1-1 
Curve Evaluation by Least Squares ? 1 
Miscellaneous 

Subroutine Library - - 
Calculation of 36 Polynomials 1 Year 40 
Center of Gravity of Parts 4 ol 
Drill Charts 8 2 


Load Distribution in Trucks 


Figure 8 


Engineering automotive problems solved by the IBM 705 
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Extracting Roots 





Eval 
It appears that the interation method by Newton 7 
func’ 
_ f(x) seri 
Xi = Xo - f!(x) 
is the best for extracting roots of any order. Let 
a 6 OS + Gg + decee + aX + ao 
whe! 
be the polynomial equation of the n‘® order, and thus the first de- ing 
rivative is divic 
vari 
f"(x) = annx™”* + an-1 (n-1)x *7+..... + @o. 1 
Inserting f(x) and f’(x) in the iteration formula results in 
x, = (Bz Lan%o + (n-2)an-1X0 | + «20. + @2X6 - ao 
.= ~ = 
RAao” + (n-l)an-1 Xo | + 000 + 2aoXo+ ay 
In order to obtain the form 
_ 
xX vu 
the constants in the equation for f(x) are 
an = 1 1 
an-1 a= 0 desi 
Se= -¥ afte: 
(n + 
Then, the interation formula becomes tain 
1 u ? 
X) -1|(n-1)x0 =| Tria 
n Xo 
I 
where u = radical, Xo = first approximation, x; = second approxi- poss 
mation, n = the order of required root. erre 
secc 
This formula allows the extraction of roots of any order and can sum 
i be programmed by the use of 39 instructions. Choosing n = 2, the func 
known formula for extracting square roots is obtained: sert 
} 
this 


1 
m= %(x0+ | J tain 
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Evaluation of Power Series 


The evaluation of exponential, logarithm and trigonometric 
functions can be made by means of power series. All of the power 
series in question have the same pattern of increasing terms: 





where C = first term of the series, B = factor by which the follow- 
ing term is increased, E, = factor by which the following term is 
divided. This factor can have the form A or A(A-1), where A is a 
variable, depending upon H. 

The factors for some series are given below. 





A B C E H 
e* 001 x 1 A 001 
sin hx 003 x* x A(A-1) 002 
Ina *) 003 x’ x A 002 
sin x 003 x" x A(A-1) 002 
sin™'x 003 x’ x A(A-1) 002 
no 


The number of terms to be used for one value depends upon the 
desired accuracy. If the answer must be accurate to seven digits 
after the decimal point, then n terms have to be used so that the 
(n+ 1)th term does not add any significant digits to the answer ob- 
tained by the summation of n terms, 


Trial and Error Procedure 


It appears that there is one general method which covers all 
possible trial and error procedures. Examples in which trial and 
error procedures are involved are polynomials of higher than the 
second order, transcendental functions, and so forth. Let us as- 
sume that the number . 1638 is the right answer of a transcendental 
function. The awkward way to approach this number would be to in- 
Sert the initial value of .0001 and increase this value by .0001. In 
this way, the computer had to execute 1,638 loops in order to ob- 
tain the right answer. The other possibility is the following: 

Insert .1000 and increase by .1000. This new value . 2000 is 
too high and has to be decreased. This can be done by shifting the 
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decimal point to the left and subtracting. Therefore, the new value 
of .1900 is obtained. This procedure has to be repeated until the 
result is found. 


. 1900 
. 1800 
. 1700 
. 1600 
- 1610 
. 1620 
. 1630 
. 1640 
. 1639 
- 1638 


It appears that this method requires only twelve loops instead of | 


1,638 loops. 


Differential and Integral Calculus 
Every differentiation or integration consists of two parts: 


1. The mathematical solution, 
2. The numerical evaluation of the solution. 


A digital computer cannot help in part 1 but it can help in part 2. 
The numerical treatment can be broken down into equations or 
series by programming for a digital computer. 

In the last four examples of subroutines, it has been shown that 
many mathematical operations, in addition to the four arithmetic 
operations, can be performed successfully by the automatic com- 
puters. 


SUMMARY 


It has been shown that the increased complexity of engineering 
problems requires new, powerful labor-saving devices in order to 
reduce the routine efforts of highly skilled technical personnel and 
in order to make possible the solution of many problems not effi- 
ciently soluble at the present time. 

Based on experience gained in actual problems, it is believed 
that the new large-scale electronic computers will meet these 
requirements. 
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Tire Vulcanization—An Application of 
Unsteady Heat Transfer 


By Samuel R. Schaffer 


United States Rubber Co. 
Detroit 


1. INTRODUCTION 


A, Vulcanization 


The change which took place when Goodyear accidently caused 
rubber to be combined with sulfur is described in today’s terms as 
a cross-linking reaction, in which the long chain rubber molecules 
are tied together by means of sulfur bonds to form a 3-dimensional 
network. A simplified picture of cross-linked or cured rubber is 
shown in Figure 1 below. Each of the chains consists of thousands 
of repeated molecular units. Goodyear termed this process vulcan- 
ization after Vulcan, the Roman god of fire, because heat was neces- 
sary to carry out the reaction. 
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Pictorial Representation of Cross Linked Rubber 


Natural rubber and most of the synthetic rubbers used in large 
volume are still vulcanized by reacting with sulfur at elevated 
temperatures. 

The most important consequence of vulcanization is that the 
rubber is converted from a plastic to an elastic material, i.e., the 
Strong chemical forces at the sulfur bonds restrain the chains from 
Sliding past each other. Therefore, cured rubber which is deformed 
under load will retract with nearly complete recovery when the load 
is released. 
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B. Effect of State of Cure on Properties of Vulcanized Rubber 


The degree or extent to which this reaction is carried is known 
as the “state of cure.” The state of cure is given by {Rd 6, where 
the rate R depends on temperature and generally varies in time. 

As the curing reaction proceeds and the state of cure increases, 
the physical properties of the rubber compound, such as tensile 
strength, tear resistance, abrasion resistance, increase rapidly, 
approach a maximum and then recede somewhat more slowly. A 
generalized curve of properties vs. cure will look like the curve in 
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Generalized Curve of Compound Properties Vs. 


State of Cure 


Since the performance of rubber items will depend on one or 
more of these properties, it is necessary that the state of cure be 
controlled within the range that provides optimum properties. In 
commercial vulcanization of most rubber products, the temperature 
of the unit is raised by applying heat to the outer surfaces and de- 

4 pending on the conduction of heat through the rubber to raise the 
temperature at the center. The outer portions are therefore subject 
to higher temperatures than the interior regions, and the curing 

reaction proceeds faster at the outside than the inside. 
It is desirable to adjust the vuicanization time and tempera- 
ture so that in the region of least cure, the cure has proceeded 
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sufficiently far so that it is out of the undercure range, while in the 
region of maximum cure it has not reached the overcure range. 

This ratio of maximum to minimum acceptable cure might be 
3:1 or 5:1, depending on the compound. The maximum and mini- 
mum values establish the first set of cure limits that must be met. 

In the case of tires, a narrower set of cure limits is required 
for the outer surface of the tread, as the resistance to cutting and 
chipping falls off rapidly with higher cure, while tread wear falls 
off with lower cure. 

For tires a third limit is imposed by the fact that either steam 
or hot water under pressure and contained within a rubber curing 
bag is used as the source of internalheat, and that some time before 
the end of the cure cycle it is necessary to exhaust the internal 
pressure down to zero so that the press can be safely opened. As 
the bag pressure is reduced, the water vapor and air which were 
trapped in the interior of the tire during the earlier assembly op- 
erations try to expand and blow the tire apart. The third condition 
then is that the cure has proceeded far enough at the time of bag ex- 
haust so that the stocks have sufficient strength to withstand the 
internal forces due to trapped gases. 


2. ONE DIMENSIONAL TRANSIENT FLOW 


A, General Case 


We have imposed a number of limits on the state of cure and 
have indicated that state of cure = {Rd ® where the rate R depends 
on temperature. For the rubber compounds used in tires, we find 
that the rate of cure approximately doubles for a 16°F increase in 
temperature. It is necessary, then, to know the temperatures that 
exist in the tire during the cure cycle in order to set up a cycle such 
that the cure falls within the specified limits. 

The physical system we are dealing with consists of a mold, 
tire and bag, with the outer mold face heated by steam and the inner 
bag surface heated by steam or hot water. The heat flow is from 
outside toward the center and can be approximated as a problem in 
one dimensional flow. To examine the conditions for heat flow in 
one direction, consider an element of thickness x and area A. The 
heat balance for this element is: 


Heat flow in - heat flow out + heat generated = heat stored. 


The curing reaction is slightly exothermic, but for soft rubber 
compounds the heat generated is small and can be ignored. (For 
hard rubber compounds much larger amounts of sulfur are used 
and the heat generated by the reaction becomes important). 
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Application of Heat During Tire Cure One Dimensional Heat Conduction 


The rate of heat flow by conduction is given by the equation: 


where 8 = time, 
k = thermal conductivity. 


The minus sign is required because the direction of positive heat 
flow is from higher to lower temperatures. Then writing the heat 
balance as a rate equation, we get: 


-kA (2) + kA (2) = A(Ax)cP 3G 
ox X) ox x, +Ax 0 


This reduces to: 


kat a 
(1) cP ax’ 20 


which is Fourier’s equation for heat conduction in one direction. 
The term is called diffusivity and is usuallydesignatedby a. 
The solutions to this equation for many regularly shaped solids, 

such as slabs, cylinders, spheres and bricks are given in a number 

of standard texts. 


















at 





TIRE VULCANIZATION 


B. Solution of One Dimensional Transient Flow Equation for 
Infinite Slab 


For an infinite slab of thickness L, initially at uniform tempera- 
ture tj, then at 8 = 0, suddenly exposed to constant surface temper- 
ature t, on both faces, the solution is: 


nx —n2M2 og @/ 2 
e 





4 - 4. 
(2) tix,0) . t, t (tj nt t,) — >> > Sin 


n= 1,3,57~~ 


This is a rapidly converging series and for most cases only the 
first term need be considered (after a relatively short time). 
For the temperature at the midplane, x = L/2, this becomes: 


4 -To o/i? 
(3) tic,e) t+ (ti-tJ ye 


From this we can See that: 


1. The time required to reach a given temperature is propor- 
tional to L’. 

2. The time required to reach a given temperature is inversely 
proportional tow. 


C. Establishment for Cure Cycles from Temperature Data 


The problem presented in curing tires can best be appreciated if 
we consider that rubber has the lowest value of a that is reported in 
Ingersoll and Zobel’s Heat Conduction. 


Compounded rubber: a .005 Ft*/hr 


Values of a for other typical insulating materials: 


. 006 (silica aerogel) to . 010 (asbestos). 


Values of a for metals: aluminum 3.3 
steel 0.48 


This means that if a set of slabs made of these materials all initially 
at the same temperature were exposed to a new surface tempera- 
ture, the rubber slab would have the slowest temperature rise at in- 
terior points. The next figure shows some temperatures measured 
during the cure of a 9.00-20 tire. This is the size that might be 
used on a city bus. 
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Figure 5 
Temperature Rise in 9.00-20 8/10 Fleetway. 
70’ McNeil U.V, Cure Cycle. 286°F External. 
310°F Circulating Hot Water 


Note that the temperatures are plotted on a scale that expands 
exponentially rather than a linear scale. The temperature scale has 
been designed so that the vertical distances are proportional to the 
rate of cure rather than the temperature. Conversion from tem- 
perature to rate of cure is based on the fact that the temperature 
coefficient is slightly greater than 2 so that the rate doubles for 
each 16°F temperature increase. 

The area under the rate-time curves gives state of cure. Using 
temperatures measured by means of thermocouples, the state of cure 
values for the surfaces and interior are calculated and checked 
against the specified state of cure limits. 

Note that in the actual case of a tire cure the temperature at 
the surface does not increase abruptly with time, so that the ideal 
boundary conditions set up for equations (2) and (3) are not met. 
How then can these temperature measurements be used to estimate 
or calculate the temperatures for a new set of conditions? 
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3. ESTIMATE TEMPERATURES DURING CURE CYCLES 


Let us consider the system consisting of the mold, tire and 
bag. Each of these has a different thickness and conductivity. How- 
ever, we can think of the mold and bag as being replaced by slabs 
of rubber having the same thermal properties as the tire. This 
system is shown in Figure 6. The outer faces of the mold and 
bag, which are exposed to steam and hot water, respectively, 
come much closer to meeting the assumption of constant surface 
temperatures than do the tire surfaces. 
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«mye, = " " “ * DIFFUSIVITY 


d,,d, = m€0uceo *' a “THICKNESS 
CORRESPONDING TO TIRE DIFFUSIVITY, X¢ 


TOTAL THICKNESS OF EQUIVALENT HEAT TRANSFER SYSTEM 
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Replacement of Mold and Bag With Slabs Having 
Diffusivity of Tire Material 
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Consider this system as a single slab, for which we wish to 
know temperatures at various depths. The general solution when 
both surface temperatures are equal is still given by equation (2). 
Using only the first term of the series, we can simplify equation (2) 
by substituting 


t \ 7 
A =5, Ax = A sin7, Bx =(—- 1)Ax, and C - <= 
os L — 2 Li 
7x -cé 
tix,@) ts , (ti - ts) A Sin = e 


(4) . a 
, ts[1 (i : ) Axe]. ts (I- Bye " 


We can write an equation of this form for any depth in the slab in- 


cluding the tire surfaces and center. 


t, = TEMP.AT DISTANCE X FROM 
SURFACE 


t > TEMP. AT SURFACE 


L= EQUIV. INFINITE SLABS THICKNESS 
2 

TA 

— 

Lé¢ 


‘© 


StLOoPE- ~Cc 
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Figure 7 


Determination of Equivalent Infinite Slab Thickness 
From Tire Temperature Data 
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Rearranging and taking logs gives: 
(5) Log (1 - *)- Log Bx - c®. 
s 


Equation (5) shows that for constant mold steam and bag water tem- 
perature a plot of a simple function of the temperature at any thick- 
ness, x, within the tire should be linear in time. 

By differentiating equation (4) with respect to time, we get: 


dt 
(6) —— = c(tg - tx). 


The temperature difference ts - tx,g can be thought of as the 
driving force tending to change the temperature at depth x; at any 
time 8 the rate of change of temperature at x is proportional to this 
driving force. This equation for temperature change in unsteady 
state conduction is analogous to the equation for heat flow in steady 
state conduction, which states that the rate of heat flow is propor- 
tional to the driving force of temperature difference. 

With temperature measurements available at various depths ina 
tire under one set of boundary conditions, equations (5) and (6) can 
be used to estimate new temperature curves and new state of cure 
values corresponding to new boundary conditions. Further, since 


and a is known or readily determined, we can find a value for L, 
which represents the equivalent thickness of the heat transfer 
system. 

In this way, from one set of experimental temperature values, 
we can predict the results under many other conditions, including 
new boundary temperatures, new initial temperatures,and changes 
in the equivalent thickness L. The final estimate is usually con- 
firmed by additional direct temperature measurements using ther- 
mocouples embedded in the tire. 
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Some Non-Parametric Methods of Testing the 
Validity of Assuming that a Sample Comes 
from a Given Population 


By Leonard G. Jonnson, Research Staff 
General Motors Corporation, Detroit 


INTRODUCTION 


Of all questions which arise in sampling statistics, one of the 
most important is whether or not it is reasonable to assume that a 
given sample comes from some definite population. Marked devia- 
tions of a set of one or more sample points from a population line 
usually create a suspicion as to whether the particular sample at 
hand really is a member of the family of samples which have a 
reasonable chance of being randomly picked from the population. 

A convenient method of determining whether or not a sample has 
a reasonable chance of coming out of a given population is to choose 
some meaningful “statistic,” i.e., some function of one or more of 
the given sample points, and then determine where the numerical 
value of the statistic falls with respect to the complete range of 
values of the statistic for all possible samples of the same size taken 
out of the population. This entails the determination of the distribu- 
tion function of the particular statistic chosen as a criterion. 


Example I - A Test Based on the Two Lowest Values 


Let us assume that a fatigue test has been run on some ball 
bearings, and we observe that there is a wide separation between the 
first and second failures, and that the succeeding failures are 
bunched close together, as shown in Figure 1. 

Let us further suppose that the line AB is claimed to be the pop- 
ulation line for this sample. We then raise the question as to the 
validity of assuming that such a sample could very readily come 
from the population represented by line AB. The wide separation 
between the first and second failures suggests the use of a statistic 
based on these two failures alone as a criterion. Let us therefore 
choose the non-parametric statistic given by the probability that the 
first failure lies below rank r;, and that the second failure lies 
above rank ry». 
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z + dz is the product, (>) g2(z)dz. But, the possible positions for z 


are anywhere from r, to 1. Hence, the probability that the first fail- 
ure lies below rank r,, and that the second failure lies above rank 
r, is 


1 


P, = (=)g2 (z)ae.. 


2 
It has been shown that for a sample of size n 
g2(z) = n(n-1) z (1-z)""? . 

(See Industrial Mathematics, Vol. 2, p. 3, Formula (A). ) 


1 


Py -| (=)n(n-1) 2 (1-2) *dz = nat (1-z)""* dz 


rs 2 


1 
n-1 
aa) = nr,(1-r,)""* ? 
Te 


0-0, 


n-1l 


which is the desired statistic to be used as a criterion. 
Suppose now, for the sample in Figure 1, that n = 50, r, = .011, 


and r, = .074. Then, for such a sample, P = 50 (.011) (1-.074)* 
= .0127. 


This tells us that, for a sample of 50 failures, the probability that 
the first failure lies below rank .011, and that the second failure 
lies above rank .074 is .0127. That is, there are about 127 chances 
in 10,000 that a sample of 50 would have the first failure below 
rank .011 and the second failure above rank .074. One might decide 
from this that the sample does not come from the population repre- 
sented by the line AB. But, how often would such a decision be 
wrong? This is a more fundamental question which should be an- 
swered before any decision is made at all. 

We can not say that the probability that this decision is wrong 
is .0127, because of the fact that even for the case in which the first 
and second points fall exactly on the line AB at median ranks, i.e., 
when r, = .0138 and r, = .0336, the value of the statistic is P, - 
50 (.0138) (1-.0336)** = .1293, and we know that there are much 
more than 1293 chances in 10,000 of being wrong when we reject the 
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hypothesis that such a “perfect sample” comes from the population In 
represented by the line AB. It therefore becomes necessary to | 
determine the distribution function of the statistic P, = nr,(1-r,)""' integ: 


for all possible samples of size n coming out of the given population, 
To do this, we proceed as follows: 
The probability density function of r2 is Thus. 


f(r.) = n(n-1)r, (1-r,)""? . 





wher 
and v 
Let y = (1-r,)""* , or rmzil ( - a. T 
2-n \ » 
n- 1 
® - ab x, 
n- n-2 
y™ rl 


io Dor 
\f{(r2)dr2| = jn(n- 1) (1-y’ *) _ ro =n(l-y~ *)dy Ther 





- hic! 
Thus, the probability density function ofthe quantity y = (1-r,)""* is oat 
F 
f(y) = n(1 - yt) 
We can now determine the distribution function of the product 
q = ri(l-r,) ' = xy as follows: (x represents r,, and y represents 
(1-r,)""*) Henc 
p32 
Prob. (a < y < a+ da) = als - ai) aa 
an r r 
- m ofe a . 
Prob. | (x <ri) c (a< y < a+ da)| = — = n-1 =. i This 
2 l-y l-a 
7 0021 
r do ae = 
Prob. fe < %) < (ac yeas da) = 1) : 
L a 1- a ee Pp 
By multiplication: or, P 
T 


prob. | (x %)- (ac ycasaa) = — 


da 
= Ndo * ° 


Another way of writing this last probability is 


Prob. Cc <qo)-(a<s y<a+ da)| = Ndo aa 








SOME NON-PARAMETRIC METHODS OF TESTING 61 


In order to determine the total probability that xy < qo, we must 


integrate nqo fs between the lower and upper limits of a. 


ae 


da a. 
Thus, Prob. ,.(xy < qo) = J Ngo — = nqo Loge|—} , 
PP . a ai 
“= 
where the subscript pp denotes the principal part of the probability, 
and where a; and a, are respectively the lower and upper limits of 
a. These limits are determined from the fact that 


-_. a . | 


a (1 - a i) 


Therefore, in the principal part, a may range over all values for 
1 


which the inequality a (sai) = qo holds. Thus, it can be seen that 


the range of integration on a is a function of qo. 
For the particular case which we are considering, we have 


*  oi27 


n = 50, and qo = 50 50 


= 000254 . 


Hence, the range of integration is determined by the inequality 


JS 
a (1 - 2 ) > .000254. 


This inequality is satisfied by all values of a in the interval 
.002157 < a < .987474. Thus, for this case, a; = .002157 and 
a, = .987474. 


. J 987474 
”. Prob. pp] rs(1-r2)" * 000255 | = 50(.000254) Log, (287474 ) 


p ~ ‘987474 - 
or, Prob. pp (Py < .0127) = .0127 Log, ao} .O778 . 


To this principal part we must add the complementary part of 
the probability, which is equal to 


a 
al 


J f(y)dy + J ty)ay = 1+ F(a,) - F(a), 
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y - ey 
where F(y) -[ n (1-y""*) dy = ny -(n-1)y"" . 


Hence, the complementary part of the probability is 
Prob. cp (P, < .0127) = 1 + F (.002157) - F (.987474) 
30 
= 1 + 50(.002157) - 49(.002157)* 


ao 
- 50(.987474) + 49(.987474)* 
.0147. 


Prob. (P, < .0127) = .0778 + .0147 = .0925. 


This last number represents the probability of our being wrong when 
we reject the hypothesis that the given sample comes from the popu- 
lation represented by the line AB on the basis that the value of P, is 
as small as .0127. 


Example II - Another Test Based on the Two Lowest Values 


Another statistic which might be used as acriterion and which 
involves only the two lowest values is the probability of obtaining a 
separation in rank between the first and second failures which is at 
least as greatash=ry, - ri. 

The formula for this probability is obtained as follows: 


leeuacriy 1 PAL/RE—pet—WO FAILURES ote ra 
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Figure 3 


The probability that the second failure is located between ranks 
z and z + dz is g,(z)dz, where g,(z) is the probability density func- 
tion of the rank of the second failure. 

The probability that the first failure lies below rank (z-h), if the 


, . Zh 
second failure lies between ranks z and z + dz is — ; 


Therefore, the joint probability that the first failure lies at leasth 
rank units below the second, and that the second is located between 
the ranks z and z + dz is the product (? -h 


) es (z)dz. 
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The possible positions for z are anywhere from hto 1. Hence, 
the probability of a separation of at least hin rank between the first 


and second failures is 


1 1 


Qn { (7-7) 2 (z)dz f 





(2:4). n(n-1) z (1-z)""? dz 


h h 
1 
nin-0f (z-h)(1-z)""? dz. 
h 

Integrating by parts: u=z-h dv = (1-z)""? az 

1- n-1 

du dz Vv _(i-z) 

n-1 





( ; | 
(z-h) (1-z)""' 1 
Q, : n(n-1) < - n- 1 |. ' Arf (1-2)""! az ¢ 


1 

n 

-ain-n oat. Eck (1-h)", which 
n-1l n h 


is the desired formula. 

For the special case considered previously, we had r,; = .011 and 
r, = .074 with n = 50. Hence, h=r, - r, = .063, 
and Q, = (1-h)"= (1-.063)*° = (.937)* = .0386. Thus, the value of 
the statistic is .0386, which represents the probability of having at 
least .063 units of separation in rank between the first and second 
failures. Now if we reject the hypothesis that the sample comes 
from the population represented by line AB of Figure 1 on the basis 
of Q, =.0386, what is the probability that we are wrong? We an- 
swer this question as follows: 

The probability density function of h is 


f(h) n(i-h)"™* . 


1 
put y (1-h)® , or h l-y" 
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| mt os Now 
sa | £(h)ah | = ny as (- 2 ay) = dy. 
*, The probability density function of the quantity y = (1-h)" is (See 
f(y) =1 (o< y <1) 
*, The cumulative distribution function of y is 
y 
F(y) f f(y)dy = y 
Oo 
Prob. (y < Qn) = Qn. 
Thus, Prob. (y < .0386) = .0386, and we see that with this crite- k = 
rion the probability of a wrong decision is the same as the numerical 
value of the statistic. Since the probability of a wrong decision is 
less with this criterion, it follows that this criterion is superior to If, 
the first one considered. This, of course, need not be true for all that 
samples, just because it happens to be true in this particular case. wha 
we 


Example III - A Test Based on the First Three Failures = 


As a third example ofa criterion we shall consider the proba- 

bility that the first failure lies below rank r,, and that the second 

and third failures are confined to a k - band whose upper boundary 

is at or above rank z. Let 
We derive the formula for this probability as follows: 


ie 
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Prob. |(3rd Failure between z and z + dz) and(exactly 1 failure| 
below r,) and (exactly 1 failure between (z-k) and z) J 


= g3(z)dz - ar ° * = arik 6, (2) dz s is 
Z 


where g;(z) is the probability density function of the rank of the third 


ark gs (z) dz over the range z =z toz=! 
Zz 





failure. By integrating 


we obtain the desired probability as 
i 


) 
Ry = af TEE es (2) da : 
Zz 
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Now, n(n-1) (n-2) 2? 


2 
(See Industrial Mathematics, Vol. 2, p. 3, Formula (A).) 


gs(z) = (1-z)""* 





R _ Bevken(n-1)(n-2) { (1-z)°"3 
a * 5 ; 





_ _2r,kn(n-1) (n-2) (aay) 
4 2 ~  n-2 


= rykn(n-1) (1-z)"~? . 
For the sample plotted in Figure 1 we have n= 50, r; = .011, 


ite- = k = .013, z = .087. 

ical 

_— .. Ry = (.011)(.013)(50)(49)(.913) = .0044. 
r to 


If, now, on the basis of this value of R , we reject the hypothesis 
that this sample comes from the population represented by line AB, 








e. what is the probability that we are wrong? To answer this question 
we must determine the distribution function of the product 
§= r,(1-z)""?. We do this as follows: 
The probability density function of z is 
ba- 
‘ond f(z) “ n(n-1) (n-2) z? (1-z)""3 
lary 2 . 
Let y = (1-z)""?, or z = 1-y™? 
3-1 
n-2 
a 
dz = a dy 
n- 3 a-2 
1 = atthe 
- It(z)d2| - | n(n-1)(n- 2) (1. yma)’ . yhme (-» ~ ay) 
2 n-2 
ure| : one) (1-y#7)* ay 
4 
Thus, the probability density function of the quantity y = (1-z)""? 
is 
' e git 12 
f(y) = n(n-1) (1-y**) : 
hird 2 
z=1 We can now determine the distribution function of the product 


§ = r,(1-z)""* = xy as follows: 


(x represents r,;, and y represents (1-z)"" 2) 
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- as We 
earner aia)" da 
Prob. (x < m)C(a<yé« a + da)| il be ma eet 
l-a™? 


Prob. I(x < §s)c @« y+] ca 
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- 2—\? Th 
-.Prob.| (x < ). (a< y < a+ da)| - a(n-1) (1-2 ny) (1- ‘ai nea) ' 
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: mnt) £8. (1087) aa 








The last probability may also be put into the form 


. A F(} 
Prob. | xy < 5 6) ‘(asc y<a+ da) | = mek (42) (1-a™* )aa, 


“ 3-2 
- Prob. ,,, (xy < §,) = tactic f (2 . af?) da 
a, 


2 


' Again, the subscript pp denotes the principal part of the probability, 
while a, and a, are the limits on a between which the following in- 
equality holds: 


. n(n-1)5o fuog.( 2 ); (n-2) (a, m3 - a 4} 





Eo <1 Fi 
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sn 
In our particular example, go 
Rn 0044 
R, = -0044, or Eo= knin-1) ~ (.013)(50)(49) ~ .000138. ok 
‘. The limits on a for the principal part are to be determined from Fr 
the inequality of 


i 
a(1 - a*) > 000138 . 
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We find that this inequality is satisfied for 
.00104< a < .99335. 


Thus, Prob. pp(& < .000138) _ 50(49) (.000138) Loge (-2232 ) 





2 .00104 
+ 48(.866667 - .99986) 
= .16905 (6.86188 - 6.39307) = .0793. 
*. Prob. pp(Rn < -0044) = .0793. 


The complementary part of the probability is given by 


ay; 1 
/ f(y)dy Lf, f(y)dy = 1+ F(a,) - Fla.) , where 
0 “2 


y 


L-\? n=) . a nh 
F(y) f n(n-1) (1-y"*2) dy ain) y - n(n-2)y™? 4 (n Ay{an3) ya 
2 
0 


Hence, the complementary part of the probability is 


Prob. ep(Rn < .0044) = 1+ F(.001040) - F(.993350) 
49 590 
1 + 1225(.001040) - 2400(.001040)* + 1176(.001040)* 
49 $0 
- 1225(.993350) + 2400(.993350)* - 1176(.993350)* 


0291. 
. Prob. (Ry < .0044) = .0793 + .0291 = .1084. 


Thus, if we decide to reject the hypothesis that the sample in 
Figure 1 comes from the population represented by line AB because 
the value of Ry, is as small as .0044, the probability that our de- 
cision is wrong is .1084. 

It can be seen that although the probability criterion R, is the 
smallest of the three criteria thus far considered, it still is not as 
good a criterion as Pn or Qn, because, as we noted before, 


Prob. (P,, < .0127) = .0925 
and, Prob. (Qn < .0386) = .0386. 


For this particular sample, therefore, the criterion Q, is the best 
of the three thus far considered, since it yields the smallest prob- 
ability of a wrong decision. 
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Miscellaneous Other Criteria and Their Distribution Functions 


I. The Criterion (h; - hz) 
Let r: = rank of first failure 
r. = rank of second failure 


r,; = rank of third failure Now 
Furthermore, let h, = the difference be rank between second and 
first failures 
‘I 


= Tg - 61 
and let h. = the difference in rank between third and second failures 
= Bs - Is 


Now take the statistic 6 = h; - he asa criterion. In order to de- 
termine the distribution function of 6 we must know the probability 
density function of h, and also the probability that h: <A if 
a < h, < a+da. The probability density function of h, =r, - r2 is [ 
found as follows: 
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Figure 5 = 1 
Prob. (Third failure lies between ranks z and z + dz) = gs (z)dz, ; n{ 
where g;(z) is the probability density function of the rank of the 
third failure. x 
° Prob. (Second failure lies between z-h, and z if the third failure 
is between z and z + dz) 
2 
= Zhe (z-he) + he - yey he) (For z > h,.) 
Zz Zz Zz 
1 (For z<h,). we | 
Hence, Prob. (Second failure lies between z-h, and z, and thatthe “ 
third failure lies between z and z + dz) , “ 
2 = he (2z-ha) gs (z)dz (For z > h2) ] 
Zz 
= gs (z)dz (For z< h,). 
be 
.. Prob. (Second failure is not more than hz rank units below the |;  o 


third failure) 


es 


dz, 
the 


the 
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he 1 
-[ gs (z)dz .f ha gs (z)dz . 
h, 2 


0 


Now g3 (z) - aa- Die“ Zz" (1-z)""3 , 


2! 
he 

~. Prob. (rz - rz < he) f ne (1-z)""? dz 
0 


i 


+f hg n{n-1){n-2) nin in (2z - h,)(1-z)""* dz 


h, 
: F - (1-z)® - nz(1-z)"™* ned) s* (i-s”” | 








hy n(n-1)(n-2) . (Qe-hy)(I-2)"* gs | * 
. 2 L n-2 (n-1)(n-2) h, 
© = (Ieha)™ = nha (1-ha)™™? ~ 2 3 (1p)? 
» MD 2 (1h)? + nh, (1-h,)"™* 
=l- (1-h,)” “ 


Therefore, the probability density function of h, is 
f(ho) = n(1-h,)""*, 


In order to determine the probability thathi<w ifa< h,< a+da 
we first prove the following: 


LEMMA: The probability that the difference in rank between the 
third failure and the first failure is not more than k is 1 - (1 - k)® 
-nk(1-k)™”. 
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Figure 6 



















70 LEONARD G. JOHNSON SOM 


Prob. (Third failure lies between ranks z and z + dz) = g3(z)dz H 
where g;(z) is the probability density function of the rank of the 





third failure. E 
Prob. [(Exactly 2 failures between z-k and z) if (third failure lies 
between z and z + dz)] 
or, | 
2 
= Ls if z>k 
Zz 
= ] if | gt ges 
Prob. [(Third failure lies between z and z + dz) and (exactly 2 
failures between z-k and z) ] 
We | 
2 
= k’ Gs (2) dz if z>k 
z” 
= g3(z)dz if z<k 
.. Prob. (Difference in rank between third and first failures is not 
more than k) 
. 1 
k? 
” f gs(z)dz + [xs gs (z)dz . 
oO 
or 
Since gs (z) = wa ine) z*(1-z)""* , this becomes 
k 
Prob. (r;-ri < k) | ain TMn-2) 92 (1-2)"* dz 
0 
1 
2 
.| k* n(n-1)(n-2) (1-2)""? az 
k 2 
(1) 


+ 
n-,; n(n-1) , 0-4] 
E - (1-z)" - nz(1-z) _—— z (1-z) 


: k ate (1-2)"4| 
K 


k? (1-k)""? 


1 - (1-k)" - nk(1-k)""* - ate) k? (1-k)""? 5 nn) 


= 1 - (1-k)" - nk(1-K™ 
6. B:. D. 
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dz Having proved the lemma, we can now put k = h; + hz. 
the Hence, Prob. (rs-ri < hi + he) 
- = 1 - (1-h; - h,)"- n(h; + h,)(1-h; - h,)™”? 
or, Prob. | (r2-ri < hi) if (r3-r, = he} 
=l1- (1-hi = h,)" aa n(hi + he )(1-hi - h,)"- A 
Prob.| (hi <u) if (ag hh, < a+ da)| 
9 = 1 - (1-H -a)" - n(ju+a)(1-u-a)"™ . 
We have, therefore, 
Prob. (a < h, < a+ da) = n(l-a)""* da 
Prob. | (hi < a2+§,) clac< hp < al 
= 1 - (1-6-2a)" - (59+ 2a) (1-80-2a)""’ 
not 


. By multiplication, 


Prob. [ (hs < a+b) (a< he ¢ a+ da)| 
= n(1-a)""* lt-(1-80-2a)™ ~ n(bo + 2a) (1-6,,-2a)""* | da 


4 


or Prob. [ (ms -h, < bo ): (asx hy c a+ da)| 
= n(1-a)"' [1 - (1-5 0-2a)® - n(§ 0+ 2a)(1-5o - 2a)" da 


Prob. pp\fu ~h2 < Bnd 


ae 


AP 
f n(1-a)"™? k - (1-8, - 2a)" - n(bo+ 2a)(1-8 9~2a)""*]a 


aj 


a2 
(1) (1-a,)"-(1-a,)"4 nin-1)f (1-a) ""*(1-6,-2a)”" da 
a, 


i2 
- of (1-a) ""*(1-5,-2a) da 


- 2 
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where a; and a, are determined from the inequalities | o*. Pi 
0 < 1 - (1-8, - 2a - n(d, + 2a)(1- 6,- 2a)" '* < 1 


For 0 <6, < landneven, these inequalities are satisfied in the 
interval 


L 1- bo ~~ 
o<a<=——. we hi 
2 
Hence, in case 0< §, <1 andn is even, we have 
ai =O 
ae = 1- bo . Henc 
2 
For the sample given in Figure 1,. we have 
h, = .074 - .011 = .063 ap 
h, = .087 - .074 = .013 
n = §0 T 
The criterion h: - hz becomes .063 - .013 = .050. -050. 


We now might ask, How significantly large is the value of hi - h, } Figu 


= .050? In other words, what are the chances that h: - hz = .050 is ae 
ever exceeded ? To answer this question we go to formula (1), and omens 
obtain (taking a, = 0 and a, = .475). ~~ 
sam] 
Prob. .,,(hi-h2 < .050) = 1 -(1-.475)" Bae 
n° 475 + 475 0 7T 
+ 2450 | (1-a)**(.95-2a)™ da - 2500 (1-a)*°(.95-2a) da. | 
0 0 

L 

The integrals in this expression are best evaluated by some ac- 

curate method of numerical integration, such as Weddle’s rule. We 

find that 

then 
“678 distr 
f (1-a)**(.95 - 2a)*® da = .000525924 prob 


0 


475 


f (1-a)**(.95 - 2a)*° da = .000492878 


0 
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. Prob. pp (hi-hz < .050) = 1 - (.525)° + 2450 (.000492878) 
- 2500 (.000525924) 
= 1- .1073 
= .8927 . 


In order to determine the complementary part of the probability, 
we have 


. 
F(y) = { n(1-y)""* dy = 1- (1-y)". 


0 


Hence, Prob. ...(hi-h, < .050) = 1+ F(o) - F(.475) 
= l+o- [1 - (.525)" | 
= .0000 + 


.. Prob. (hi-h, < .050) = .8927, or Prob. (hi-h, >.050) = .1073. 


cp 


This last number tells us the significance of the value h, - hz = 
.050. We see that if we reject the hypothesis that the sample of 
Figure 1 comes from the population represented by line AB, the 
probability that we are wrong is .1073, since it is theoretically pos- 
sible for 10.73 per cent of all samples of size 50 drawn from the 
population to have a value h; - h, > .050. Hence, for this particular 
sample, the criterion 6 = h; - hg is also inferior to the quantity Q, 
of Example II. 


I. The Criterion A = D_ (l< js n-1) 
Tj +1 
Let 5 = rank of the jt® failure 
r+, = rank of the (j + 1)th failure, 


J 


then consider the statistic defined by the ratio A = i. . what is its 
rj +1 

distribution function? To answer this equation we must know the 

probability density function of r;,,, and the probability that 


ry¥< Ao if a< Ti < a+da. 
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It has been shown that the probability density function of Ti+. is 


g(z) = 


where z represents Pj+s 





j! (n-j-1)! . 


ni 


j (1-z) 7", 


(See Industrial Mathematics Vol. 2, p. 3, Formula (A). ) 


Furthermore, it can be seen from Figure 7 that 


Prob. | (Exactly j failures below rank.) if 


Thus 


f “Es r i FAILURE 
y } ” RAS ——m 2 4+ 
£4 ae gr 2 puver 
z s | i a 
A, 2 z+de 
Figure 7 
Prob. (as r,, < a+da = a ee a’ (1-a) "F"" da 
rob. je. 7 j!(n-j- )! 
r< 
and Prob. lin) < A,)Cc lag rjti < a+ aa) » 4° 


o 


(j + 1)t2 failure lies between z and z + az)| = — ; 





J 
Prob. ey ad.) Cc (a €Tj+i | a+ da)| =Ao. 


By multiplication: Prob. [(r, < ado)- (ac Titi € at da) | 


~ j! (n-j-1)! 


(25° ho) ‘(ag Pj+i Ss a aa 
Py +1 








n! 


Aj a? (1-a) "2 da 


Aga (1-a)™* da 


sAZ . 
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— r 
Hence, the statistic A = —/ has a cumulative distribution 


"jas 

function F(A) =A, i.e., a probability density function f(A) = jx” *. 
This criterion \ turns out to be very useful in case two succes- 

sive upper extreme values are widely separated. Let us say, for 

example, that in a sample of 51 specimens the 50th failure is at 

rank .920, and the 51st failure is at rank .980 of a given population. 

Then, taking j = 50, 


a=3 i oe = 9387755 
rj+i Is . 





and Prob. (A < .9387755) = (.9387755)*° = .0425. 


This indicates that such a sample would not readily come from 
the population in question, since a value of A as lowas .9387755 is to 
be expected only in about four samples out of a hundred. 
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Analytical and Experimental Studies of 
Vehicle Dynamic Behavior * 


By John Arthur Carlson 

Mechanical Engineering Research Department 
Armour Research Foundation of Illinois 

Institute of Technology 

Chicago 


SUMMARY 


Progress in the development of analysis in the field of vehicle 
dynamics has been made possible by modern computational tech- 
niques. In this study, an appraisal of analysis is made by comparing 
the results of analysis with actual vehicle tests, both experiment and 
theory being included in one investigation. The general conclusion 
drawn from such an appraisal is that relatively simple systems can 
be used to represent a vehicle in motion over a variety of courses 
to give good quantitative results; and further, the effects of varying 
vehicle properties (spring rates, damping characteristics, masses) 
can be studied by using these simple systems, the results obtained 
being of sufficient accuracy to be used for design purposes. Con- 
clusions reached on the effects of some automobile properties (sus- 
pension bump stops, suspension damping, and tire damping) are 
summarized. 


INTRODUCTION 


As the ability of ground vehicles to attain and sustain higher 
speeds has increased in recent years, the dynamic behaviour of such 
vehicles has become more critical, and the problems of passenger 
comfort and safety, and vehicle stability and control have become 
more complex and difficult. Basically, the problem is to design a 
vehicle to satisfy the requirements of the personnel or material 
transported. The goal is to obtain the optimum solution to this 
problem in the most economical way. Analytical design techniques 
which may be applied to the vehicle to predict its behaviour are 
required. 

In the work described here, the problem of obtaining the dynamic 
response of a vehicle was studied with the above goal in mind. The 
literature on this subject reveals that many analyses have been 
made from which general conclusions regarding the dynamic 


*Based on a thesis submitted in May, 1955, in partial fulfillment of the re- 
quirements for the degree of Doctor of Philosophy at the California Institute 
of Technology. 
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behaviour of vehicles have been drawn.’’”’” Practically all of the (3) 


analyses are limited to linear systems, and in only a very few cases 

is the response of a system computed. Most papers involve only, (: 
discussion of the qualitative effects of various terms which appear , 
in the equation of motion, the discussion sometimes being based op 

only the one degree-of-freedom system concept. In practice, an at- 

tempt is made to apply the meager analytical conclusions available. (4) 
however, most designs evolve through testing. With the advent of 

modern computational techniques, it is now possible to reappraise 

the problem of analyzing a vehicle and predicting its behaviour.*’® 


VEHICLE TESTING AND COMPUTER DEVELOPMENTS 


The development of modern instrumentation has made possible 
accurate and extensive determination of vehicle dynamic behaviour 
from testing.°’’’*’® While much information can be gained ex- 
perimentally, it may be seen that purely experimental design studies 
in which vehicle properties are to be varied must be limited within (5) 
economic bounds. These limited studies give an incomplete picture 
and sometimes lead to inconclusive results. 

Another approach to the problem of vehicle design is to analyze 
the vehicle as a mechanical system of many degrees of freedom 
which represents the motions of the vehicle components, such as 
wheels, frame, body, engine, and passengers. The troubles which | 
face a theoretical analysis can be seen by considering an example. 
A system which can be used to represent an automobile and its sus- 
pension approximately for small vertical motions is shown in Fig- 
ure 1. Many simplifying assumptions were necessary to arrive at (6) 
this system, such as the assumption that the engine and the body and | 
its contents are rigidly connected to the frame which is rigid and 
that the tire contacts the ground at a point directly below the wheel. 
The equations of motion of this system are as follows: ' 


( Aaa oe 4 

mi Xi -Gi [Xs + ¢ é+a6-in) - Ki (x + d 6,24-x,) e Ke (4 ‘a rw 

(1) + Ks(xi-yi) = 0; 
, d “— d i 
M2 X2 -G2 (i - d § sai ) - K, 4 -  6.a-x: ) 
Ky ( Xi —*: ) 
2 — (@- ——= } + Ke(x.-y2) = 0; 

( ) + q q 4 27Y2 


m I, \- ms; i \.. Eee e 
(3) «(Boal fib) ml fone 
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(3) - = 6- =) + Ky(xs-ys) = 0; 
(= -2)is + (= +a % -G (i . 9 -b$-i) -K, (xs -£6-bp-x,) 
(4) + ie (0- aS) + Ke (x4-ys) = 0; 


m, Xs G1 (is +4 + c¥ ~ 64a-x x1) +G, (x -¢ - — §+a- i) +Gs (xs +£ 8- bé-is) 
; e+ O « d 
+ Gs (cs - < §-pd-ix) +Ki (xs +2 6+ap-x1) 
( d 
+ K. Xs5 - 9 O+ ad-X2 
alee S esate) me fhe) = 0 
. —— ° “fe ae are 
lig + aG (i, 9 osap-i:) +aG2 (i, -4 bsg) 
- 06s (is + $ 6-borix) -bGa (5 - £6-bp-i.) 
+ aK, (x. + + + a-x) +aK, (xs - $ @+ ap-x2) 
(6) - bKs (xs +5 6-bh-xs) -bK, (x, : € 4-bd-x,) 
- ; a = d : oo 
1,6 + oG: (i. + 4 § .ag-i:) - 9 G2 (x, - $ j.ag-k) 


+26, (i hats 


+ 5 Ks (x. + § += @- -vp-xs) -$ -= Kz (xs - £ 9-b$-x.) 


os (BB) ae (BEB) «0 
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Obtaining solutions to these equations by conventional analytical 
methods would not be feasible unless many more simplifying as- 
sumptions were made.*® As simplifying assumptions are made the 
solutions may become increasingly unrealistic and the information 
obtained becomes more uncertain. The next step is to turn to the 
use of an electrical computer; but this is not as simple as it sounds, 
Sclving the above equations, which by no means represent the vehicle 
completely, with the addition of complications such as bump stops 
and non-linear shock absorbers, would tax the capacity of even the 
largest computers available. Even if solutions to these equations 
were obtained, one still would not be certain of how well the actual 
vehicle would have been represented. These computers are very 
costly and require specialized personnel to operate them. That both 
of these factors tend to hinder their ready use by engineers is a 
point which may be overlooked by computer designers. Besides 


REAR 














Figure 1 


Seven Degree-Of-Freedom System for Representing Small 
Vertical Motions of an Automobile: m, - Chassis and Body; 
m;, M, - Front Wheels; m, - Rear Axle; K;, K2- Front Sus- 
pension Springs; K,, K,- Rear Suspension Springs; Ks, K,g,Kz7, 
Ks, - Tires; G - Shock Absorbers; Ky, K,,. - Roll Stabilizers. 
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this, the given data may not be known with sufficient accuracy, or 
extreme accuracy may not be necessary for design purposes, to 
justify an elaborate computer setup. 

In an effort to resolve the above difficulties in testing and analy- 
sis, a research program was set up, the main object of which was 
to conduct an appraisal of analysis by comparing the results of 
analysis and actual vehicle tests. It was hoped that simple mechan- 
ical systems (simplified over that represented in Figure 1) could be 
used to represent the vehicle and compute its response to a given 
excitation. A previous work along this line“ was very limited in 
scope, and the analysis was performed independently of the vehicle 
tests. The work described here is the first known to the author 
which included performing the vehicle tests as well as the analysis. 
Including both experiment and theory in one investigation allowed a 
complete study of all methods and procedures used. This is a neces- 
sity if the evaluation of the analysis and vehicle tests is to be valid. 


Accelerometers 


f 
- ee 


Hinge bolt 


Figure 2 


Installation of Accelerometers at the Right Front Wheel. 
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Hinge bolt 
ww 





Figure 3 


Installation of Displacement Pickup at the Left Front Wheel. 


VEHICLE TESTS 

A standard production automobile was instrumented so that the 
vertical accelerations of the frame and axle at each wheel, and the 
relative vertical displacement between each wheel and the frame 
could be measured. Using acceleration as a measure assured an 
adequate description of the response, since velocities and displace- 
ments can be obtained accurately from the acceleration by integra- 
tion. Further, the significance of any higher derivatives in vehicle 
rideability is not well established.“ The installation of accelero- 
meters atthe right front wheel is shown in Figure 2. Figure 3 shows 
the displacement pickup installed at the left front wheel. Accelero- 
meters and the displacement pickup installed at the right rear wheel 
are shown in Figures 4 and 5. Speed and distance travelled were 
measured by means of a fifth wheel. The arrangement or recording 
equipment within the test car is shown in Figure 6.* 

Courses were laid out over which the vehicle was driven ina 
straight line at constant speeds. One of these courses contained a 
large dip which excited low frequency symmetric motion of the body 
(bounce and pitch motions). A profile of this course is shown in 


*For a complete description of the instrumentation, see the thesis referred 
to in the note at the beginning of this article. 
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Accelerometer 


Displacement 





Figure 4 


Displacement Pickup and Accelerometer Mounted on the Frame 
at the Right Rear Wheel. 





Figure 5 


Displacement Pickup and Axle Mounted Accelerometer at the 
Right Rear Wheel. 


83 





84 


Level (feet) 


0 


wi 


.0 





JOHN ARTHUR CARLSON 


j } Recorder 


Amplifiers 





Displacement Pickup 
Control Box 





Figure 6 


Arrangement of Recording Equipment in the Test Car. 
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Figure 7 


Profile of Test Dip. 
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Figure 7. Higher frequency motions of the wheels (wheel bounce) 
were excited by driving the car over short bumps of wood which 
were placed on a smooth pavement. These are shown in Figure 8 in 
position to excite right and left wheels at the same time (symmetric 
motion). By removing one of the bumps the wheels on only one side 
of the car could be excited. 

Records of the accelerations of the frame and wheels, and of the 
relative displacements of the wheels with respect to the frame, were 
made while driving the test car over these courses. Figure 9 shows 
a typical record made in traveling over the test dip. Typical vehicle 
test records obtained while going over one of the wooden bumps are 
shown in Figure 10. 

A four degree-of-freedom system which could be used to repre- 
sent the motion of the vehicle in symmetric motion is shown in 
Figure 11, together with the differential equations of motion which 
describe it. 

For the vehicle used, the front to rear mass coupling term Mi, 
(see Figure 11) was found to be small compared with the other mass 
terms Mi; and M22, which is typical of most modern American cars. 
If this term is neglected in the analysis the system shown in Figure 
11 reduces to two two degree-of-freedom systems. This two degree- 
of-freedom system, shown in Figure 12 together with its differential 
equations of motion, was used to obtain the analytical results for 
bounce and pitch motions. 

An examination of Figure 10 reveals that wheel bounce of each 
of the front wheels is practically independent while the motion of 
left and right rear wheels is strongly coupled. A system which was 
used to compute the motion of the rear axle in higher frequency 
asymmetric motion is shown in Figure 13 with the corresponding 
differential equations of motion.* It is assumed that this motion of 
the rear axle is independent of the higher frequency motion of the 
front wheels. 

To obtain solutions to the differential equations of motion of the 
above systems for arbitrary input functions in a reasonable time, 
computers are a necessity, especially if the equations are non- 
linear. The computed results presented in the next section were 
obtained by using a small electric analog computer, the principles 
of operation of which are presented in the literature.*’“ Operation 
of this type of computer is based on the fact that the mathematical 
equations describing the reaction of electrical circuits containing 
inductors, resistors, and condensers, are the same as those for 
analogous mechanical systems. This type of computer is very suit- 
able for studying the dynamic behaviour of mechanical systems. 
Other kinds of computers which can be used are the electronic dif- 
ferential analyzer (sometimes called “electronic analog com- 
puter”),*® the mechanical differential analyzer,*® and the digital 
computer.” 


*For derivations of the equations in this section, see the thesis. 
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Figure 8 


Side View of Wooden Bumps Placed Symmetrically. 
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Figure 9 


Typical Vehicle Test Record Obtained While Travelling Over 
the Test Dip at 25 mph. 
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Figure 10 


Vehicle Records Obtained in Going Over a Wooden Half-Sine 
Bump With the Left Wheels Only, at 15 mph. 
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Four Degree-Of-Freedom System for Representing an Auto- 
mobile in Symmetric Motion. 
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Two Degree-Of-Freedom System for Representing an Auto- 
mobile in Symmetric Motion. 
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Four Degree-Of-Freedom System for Representing an Auto- 
mobile in Symmetric Motion. 
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Figure 12 


Two Degree-Of-Freedom System for Representing an Auto- 
mobile in Symmetric Motion. 
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Figure 13 


Four Degree-Of-Freedom System for Representing Motion of 

the Rear Axle: m, - Chassis and Body; m;- Rear Axle; Ks, 

K, - Rear Suspension Springs; K,, K, - Tires; G3, G,- Shock 
Absorbers; Kio - Roll Stabilizer. 


2 K, K, 
q--iP tiptoe, +(G, + Mic, +x 3+ St, +x 2tay)= 0. 
a *p 
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VEHICLE PROPERTIES 


In order to compute solutions of the equations of motion it is 
necessary to know the physical characteristics of the vehicle, 
masses, moments of inertia, spring rates, and damping. For the 
vehicle used in the tests described here, these properties were ob- 
tained from data supplied by the manufacturer or from laboratory 
measurements. 

The front and rear suspension spring force-deflection curves 
were measured and are as shown in Figures 14 and 15, respectively. 
In the analysis, these were approximated by the dashed lines shown. 
Approximating the shock absorber force-velocity curves by straight 
lines and adding the average static friction of the suspension ob- 
tained from the suspension spring data, yields the damping charac- 
teristics shown in Figures 16 and 17. The force-deflection curves 
for the tire were as shown in Figure 18. In the analysis these curves 
were approximated by the dashed lines. 


RESULTS 


Low Frequency Motion of the Sprung Mass 


The results of the vehicle tests and of the analysis are presented 
in the form of graphical comparisons of the responses of the vehicle 
obtained by the two methods. 

Linear systems were used first in the computer studies in order 
to determine how much could be learned under this very radical 
simplifying assumption and also to provide a basis for comparison 
with non-linear systems subsequently studied. The assumption that 
the front-to-rear mass coupling term could be neglected when con- 
sidering symmetric motion was checked on the computer and found 
to be valid. 

Records showing the response of the vehicle in going over the 
test dip are reproduced in Figures 19 to 22. The measurements 
used in these comparisons are the relative displacements between 
the frame and the wheels, and the accelerations of the sprung mass 
m,. Low frequency motion of the sprung mass on its suspension 
predominated, and higher frequency motion of the unsprung mass as 
well as noise present in the frame, largely due to the engine, were 
neglected. The vehicle measurements plotted are an average of 
those made on the left and right hand sides, since the motion was 
not exactly symmetric. The computer solutions appear only where 
they depart from the vehicle measurements. In these records the 
Only non-linearity taken into account was that of the change in 
Spring rate occurring when the main suspension spring bottomed on 
the rubber bump stops. Tire pressure was 24 psi gage, and the 
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values of viscous damping used for each wheel in the front and rear 
were 44.0 slug/ sec. and 47.8 slug/sec., respectively. 

The agreement between the computed results and the vehicle tests 
for motion over the test dip is surprisingly good if one considers all 
the sources of error that exist between the two curves. These in- 
clude errors in the vehicle test measurements, measurements and 
calculation of vehicle properties, approximation of vehicle prop- 
erties, surveying of the course, and overall computer error. 


Effects of Changing Vehicle Properties on Low Frequency Motion 


From the above records it may be concluded that the most 
influential non-linearity of the vehicle suspension in low frequency 
motion is that of the main suspension springs caused bythe presence 
of the bump-stops. The effects of changing the characteristics of 
these stops was studied on the computer. In Figure 23, the pro- 
nounced effect of changing the rebound stop clearance on the lower 
peaks in Figure 19 can be seen. Varying the ratios of the bump- 
stop stiffnesses to the suspension spring stiffnesses affects the peak 
displacements and accelerations as shown in Figure 24, these re- 
sults again being obtained on the computer. 

Typical records obtained by using the computer are shown in 
Figure 25. These illustrate the effect of varying the main suspension 
spring rate on response motion when going over the test dip. Both 
the amplitude and frequency of the motion can be seen to increase 
with an increase in spring rate. The curves for a suspension spring 
rate of 106.2 lb/in. correspond to the non-linear computer curves 
shown in Figure 20. (Note that the direction of travel is reversed.) 

Evidently the non-linearity of the suspension damping is not 
important in primary motion of the frame since linear viscous 
damping gave satisfactory results. The values of viscous damping 
used were chosen by approximating the damping characteristics by 
straight lines in the range of relative velocities between the wheel 
and frame indicated by the slopes of displacement-time curves 
obtained in the vehicle tests. Figure 26 shows that increasing the 
amount of viscous damping lowers the amplitudes of the motion. 
This would be expected since examination of the excitation when 
going over the test dip reveals a large component with a frequency 
close to the natural frequency sprung mass on the suspension 
springs. The conclusion reached is that for low frequency excita- 
tion the damping of the sprung mass should be close to critical. 

Neglecting tire damping in the analysis appears to be justified 
by the results in the case of motion over the test dip. 
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Effect of Varying the Rebound Stop Clearance at the Front 
Wheels in Going Over the Test Dip at 20 Miles Per Hour. 
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Figure 24 
Effect of Varying Bump Stop Stiffnesses at the Front Wheels 
in Going Over the Test Dip at 25 Miles Per Hour. 
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Figure 25 


Effect of Varying the Main Suspension Spring Rate on the 
Motion at Front Wheels in Going Over Test Dip at 25 Miles 
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Wheel Bounce 


Records obtained for motions over the half-sine bumps are 
reproduced in Figures 27 to 31. Measurements chosen for compari- 
son are the wheel accelerations and the displacements of the wheels 
relative to the frame. In Figures 27 to 30 for symmetric motion, 
computed results are shown for both linear and non-linear systems. 
The values of viscous damping used in the linear cases were 44.0 
slug/sec/ wheel and 90.5 slug/sec/wheel at the front and rear, re- 
spectively. Non-linearities included in the computed results were 
the effect of tires leaving the ground and suspension damping (Fig- 
ures 16 and 17). The actual vehicle measurements plotted for 
symmetric motion are the averages of those obtained on left and 
right hand sides of the vehicle. While the agreement is not as good 
as was obtained for the dip, it is still quite good considering the 
possible errors involved. In addition to the sources of error men- 
tioned previously, here the tire was assumed to contact the ground 
at a point and the suspension bump stops were neglected in obtaining 
the analytical results. 

A typical computer record of motion of the vehicle is shown in 
Figure 32. It can be clearly seen that the coulomb friction in the 
suspension gives rise to high rates of change of acceleration of the 
frame. A passenger in the vehicle commonly refers to the quality 
of the resultant ride as “jerky.” Quantitative measurements of the 
magnitude of this effect could be made from computer records such 
as this. The damping force curve which has greater peaks in re- 
bound reflects the unsymmetrical shock absorber characteristic. 

No significant change in the primary motion over the test dip was 
found to take place with variation in the tire spring rate (tire pres- 
sure). However, varying the tire spring rate did have a strong effect 
on the peak accelerations experienced by the wheels as shown by the 
computed results in Figure 33. Higher frequencies of secondary 
motion were also observed with increasing tire spring rate,as would 
be expected. 

Figure 31 shows results obtained for motion of the rear axle 
going over the wooden bump with the left wheels only. Analytical 
results are presented for a linear system only. It may be concluded 
from the results for symmetric motion that the inclusion of the 
non-linear damping and the effect of the tire leaving the ground in 
the system used for studying the rear axle motion would yield re- 
sults showing agreement similar to that of Figure 30. 


Tire Damping 


The effect of the hysteresis of the tire, exhibited by the static 
load-deflection curves in Figure 18, was approximated in the 
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Figure 26 


Effect of Varying Viscous Suspension Damping at the Rear 
Wheels in Going Over the Test Dip at 25 Miles Per Hour. 
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analysis by adding viscous damping, G2, to the wheel. Values used 
for this damping were determined by trial and error. Computer 
results using different values were compared with vehicle test 
records until the best agreement was obtained. It was found that the 
values thus determined variedconsiderably with speed of travel and, 
consequently, with the amplitudes of the motion. 

If one looks at the tire load-deflection curve for a tire pressure 
of 24 psi gage and assumes thatthe hysteresis loop is due to coulomb 
type friction (i.e., the width of the hysteresis loop remains constant 
regardless of amplitude), the magnitude of this friction force may 
be estimated as one half of the difference in load across the loop at 
a given deflection. The static deflection of the tire was approxi- 
mately one inch; the coulomb friction force estimated from this 
portion of the curve is about + 40 pounds. The equivalent viscous 
friction “ g which will effect the same dissipation of energy as 
coulomb friction +F for sinusoidal motion with an amplitude A and 
a frequency w is given by 


g = 4F/7Aw. 


This equation may be derived by calculating the amount of energy 
dissipated by each form of damping in one cycle and equating these 
energies. For motions of the unsprung mass, the natural frequency 
may be calculated from the vehicle test records. From the acceler- 
ation records in Figure 10, the period of this motion is about 0.11 
second so that w = 27/0.11 = 57 radians/second. The amplitude of 
the relative motion between the wheels and the frame for the first 
cycle of vibration after going over the half-sine bumps may be used 
for the amplitude A. This amplitude may be read from the graphs 
of the vehicle records. Using values of A obtained from Figures 27 
to 29, and the values of F and w given above, the values of g given 
in Table 1 were calculated from the above equation. Also listed are 
the values of tire damping used on the computer. 


TABLE 1 


Comparison of Calculated and Experimental 
Values of Tire Equivalent Viscous Damping 





Speed Amplitude Equivalent Viscous Damping 
A a 
Calculated Computer 
(miles/hour) (inches) (slug/ sec) (slug/sec) 
10 0.2 54 57.5 
15 0.75 14 20.4 





20 1.4 7.7 2.6 
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Figure 31 
Motion of Rear Suspension in Going Over Half-Sine Bump With 
Left Wheels Only at 20 Miles Per Hour. 


Considering the rather crude assumptions involved, the agree- 
ment between the values calculated and those determined by trial 
and error is good enough to conclude that the tire damping can be 
approximated from measurements of the static friction. Also, the 
agreement substantiates the statement that the damping added to 
obtain agreement between the computer results and the vehicle tests 
should be attributed to the tire. 
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Figure 32 


Computer Record Showing Motion at Front Wheels, Over Half- 
Sine Bumps at 20 Miles Per Hour. 


CONC LUSIONS 


During the course of conducting the research program from 
which the above results were drawn, the author reached certain 
conclusions on the effects of some automobile properties, which 
may be of interest. 

1, The maximum accelerations experienced by the frame are 
strongly influenced by the presence of the suspension bump stops. 
Decreasing bump stop clearance and increasing the stiffness of 
these stops both have adverse effects. 

2. Neither the non-linearity of the suspension damping character- 
istics nor the unsymmetrical action of the shock absorbers have 
significant effects on low frequency (about 1 cycle/sec) motion 
of the sprung mass. 

3. Coulomb friction in the suspension system has a considerable 
effect on the accelerations of the frame caused by motion of the 
wheels at higher frequencies. 

. The incorporation of damping between the wheels and the frame 
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Figure 33 


Effect of Varying the Tire Spring Rate at the Front Wheels, 
Half-Sine Bumps, 15 Miles Per Hour. 


is wrong in principle. Ideally, the absolute motions of the sprung 
and unsprung masses should be damped independently. The prac- 
tice of damping the relative motions between these masses is a 
compromise, at best. As a result, there is a wide variation of 
opinion about what the characteristics of shock absorbers should 
be. 

. Damping in the tire appears to have the characteristics of cou- 
lomb type friction. Therefore, the effect of this damping is 
greatest for small amplitudes of wheel motion, and the effect 
decreases as the amplitudes increase. Whether or not this damp- 
ing may be neglected depends on the road excitation. 


Of greatest importance is the following general conclusion: 

Relatively simple systems can be used to represent the vehicle 
in motion over a variety of courses. These simple systems can 
be used to study the effects of varying vehicle properties. The 
results obtained are of sufficient accuracy to be employed for 
design purposes. It is possible to study the effects of different 
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properties (spring rates, damping characteristics, masses) sep- 
arately. The practically unlimited number of variations make 
studies possible which will yield a clearer picture of the phenom- 
ena taking place. 
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Improving the Validity of Transducer Measurements 
with Computing Circuits 


By S. A. Lippman 
Tire Division U. S. Rubber Company 


The meaning of the word "transducer" has undergone a change 
in recent years. Formerly transducers were devices for converting 
energy from one form to another. Now they are also instruments 
for creating variations in one physical quantity equivalent to varia- 
tions in another. The secondary variation is usually generated for 
reasons of convenience. For instance, in the study of sounds, a 
transducer called a microphone produces electrical signals which 
correspond to accoustic variations in air pressure. If we wish to 
record the character of the sounds or to analyze frequency spectra, 
we will find that the operations are easier with the electrical signals 
than with the pressure variations. In contrast to the earlier defini- 
tion, some transducers utilize sources of energy other than the 
primary variable. The goal is often to attain changes in some phys- 
ical quantity at a large power level, while drawing very little power 
from the system under observation. 

In this paper we will be concerned with the use of transducers 
for translating mechanical variables into electrical signals. Veloc- 
ity, acceleration, pressure and displacement are examples of some 
of the primary variables. In all cases the conversion to electrical 
signals requires the displacement of some sensing element. In 
piezo-electric detectors there is the distortion of a crystal; in in- 
ductive pickups there is a moving armature or coil; and in strain 
gages there is the extension of the wires and motion of the supporting 
structure. Similar motions occur in differential transformers, 
potentiometers, variable capacitors, magnetostrictive heads and 
others. The requirement of a moving mass in these transducers 
results in limitations in their use and distortions in the process of 
conversion. A somewhat general method for coping with these lim- 
itations and distortions will be the theme of this paper. However, 
before undertaking this phase of the presentation, let us consider an 
elaboration of the meaning of "transducer. "' 

Transducers are usually employed to study the behavior of a 
System when it is reacting to some kind of excitation and are also 
used to study the excitation that produces a known response. An 
accelerometer, for example, may be used to evaluate the forces 
acting on a vehicle. In this instance the output of the transducer is 
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not only determined by the variable being studied (the forces which 
act on the vehicle) and the properties of the accelerometer, but 
also by the system to which the device is attached (the vehicle). Ip 
a sense the entire system including the transducer itself may be 
treated as a more complex single transducer. Regardless of the 
terms we use, however, we must account for the influences of the 
transducer and the system being studied in order to have valid 
measurements. 

The basic idea of this technique requires knowing the differential 
equation of detection by the overall transducer. Knowledge of the 
coefficients, however, is not required. The differential equation is 
expressed in terms of the variable being investigated and the output 
signal. A circuit is then constructed that simulates this differential 
equation. Next, the circuit and the transducer are matched to each 
other in a test procedure. We now apply the output of the transducer 
.o this circuit and produce a new signal which is equivalent to the 
variable being studied. Let us illustrate the process by means ofa 
concrete example. 

The diagram of Figure 1 indicates the basic mechanism of a 
typical seismic transducer. This instrument may be used for the 
detection of displacements, velocities, or accelerations. The sens- 
ing unit shown in the diagram may be a strain gage, a slide wire, a 
reluctance bridge, a differential transformer, etc. The entire in- 
strument moves in the process of its operation. An internal mass 
M, is suspended from the case of the instrument by means of some 
flexible member. As the case moves, inertia tends to keep the 
mass at a constant position or at a constant velocity. This tendency, 
however, is compromised by the flexible member. In addition, 
there may also be a viscous drag which acts between the moving 
mass and the case. This drag is often purposefully included to min- 
imize oscillations but may be due to mechanical hysteresis of the 
flexible member or to the "radiation" of vibratory energy out of the 
transducer to the system being measured. Finally, the operation of 
the transducer also depends on the sensing unit. This device detects 
the displacement of the mass from its position of rest in the case or 
in some instances it detects a property of the relative motion. 

When we consider the quantitative relationships in the trans- 
ducer, we obtain the following differential equation: 


(1) M(X + Y¥) + KY +yn¥ 0, 
where X is the displacement of the case, Y is the displacement of 


the sprung mass (M) relative to the case, K is the spring constant of 
the flexible member, and uw is the coefficient of the viscous drag 
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Schematic Diagram of Seismic Transducer 
ent of 


F The transducer moves an amount X from the reference position. 
ant ol 


The corresponding movement of the sprung mass from its in- 
drag ternal reference is Y. The sensing device detects some prop- 
erty of the variation Y. 
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Let us first assume that the sensing unit is such that the output 
signal is proportional to Y and that the constant of proportionality is 
1. With very little loss in generality, then equation (1) relates the 
output signal to the input motion. 

In order to see how this one device can act variously as a dis- 
placement gage, velocity detector and an accelerometer, we will 
rearrange equation (1) and integrate it twice. We obtain the follow- 
ing equations. The constants of integration are omitted since we are 
only interested in variations from a reference condition. 


(2) X = -(¥Y + Yu/M + YK/M) 
(3) x = (+ Yu/M + © fyat) 
u K 
(4) K=-(¥+ & [vat + 5 f frat at) 


When the nature of the input and the coefficients are such that the 
first two terms in the brackets of equation (2) can be neglected, the 
output is almost proportional to the acceleration. For another set of 
values and functions, the first and third terms in equation (3) are 
negligible and the output is almost proportional to the velocity. Sim- 
ilarly, we can see from equation (4) that there are realms in which 
the transducer approximates the displacement. 

It may also be of interest to note the exact form of the output 
signal. If the accelerometer has been at rest at time t = 0, then 
integration of equation (1) yields: 


t 
(5) Y=- if x (t-r)e" #7/2M Sin CrdrT, 
Oo 
f 2 
where C=1/2 (a 
YM ™ 


T is a variable of integration, and the brackets indicate that X is 
taken as a function of t- 7 during integration. This reiation may also 
be expressed as the sum of the desired derivative and an error 
term. However, for reasons of brevity I have not included these 
forms. The relations of equation (5) are, in general, too complex 
for easy discernment of the details of X from the output signal. 

Now let us look briefly at the frequency response of the trans- 
ducer. This is a common way to evaluate the performance. 

For a sinusoidal input at an angular frequency w, the ratio of 
input to output is termed the transfer function T(iw) and is given by 
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(6) Tw) - ax 

wM~ Mw” 
In different ranges of performance, the transfer function can be 
almost equal to 1, creed or - Met depending on the relative values 
ofw, , M and K. Again, in each range the device acts as a dis- 


placement gage, a velocity detector or an accelerometer. In any 
frequency range the transducer is only somewhat correct in its con- 
versions and is likely to be in considerable error if the subject 
matter covers too wide a frequency band. 

The solution of restoring the desired derivative of X from the 
output of the transducer is contained in equations (2), (3) and (4). 
We have only to form the appropriate derivatives and integrals of Y 
and add these in the correct proportions. Graphical and numerical 
methods are possible, but are cumbersome as well as inaccu- 
rate. The use of an electrical analog computer at the time of meas- 
urement, on the other hand, results in a continuous flow of cor- 
rected data. The next point to be covered, then, is the nature of 
electrical analogs that might prove to be satisfactory for our pur- 
poses 

There are many approaches to the design of the electrical ana- 
logs. We might avail ourselves of the transconductance properties 
of vacuum tubes. In this way we can cause currents which are pro- 
portional to the transducer output to flow through networks. The 
resulting voltage drop is equivalent to the computed function. We 
might employ the current transfer characteristics of transistors in 
a similar manner, and we might even resort to the use of passive 
circuits for restricted frequency ranges. Finally, there is also the 
use of feed-back circuits—the method currently popular in com- 
mercial analog computers. Figures 2 and 3 contain schematic dia- 
grams for a number of these approaches and indicate some of the 
quantitative relationships. The choice of a particular type of circuit 
is not important to our discussion. In all cases the use of the analog 
computer involves the same subsequent steps. 

The constants of the transducer are presumably unknown, and in 
designing the circuits we do not usually match the values of the cir- 
cuit elements to the exact requirements. Instead, adjustments are 
provided so that a sufficient number of parameters can be varied at 
will. The scheme of matching the computer to the transducer is 
based on a test procedure which taxes the accuracy of the computer. 
The input must be a known function of time, and incorrectly com- 
puted outputs should be readily recognized. The input might be step 
functions or sharp pulses (delta functions). The response of a poorly 
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enable the tube to act as a 
capacitors are employed, we 
(CCD) to provide the operating 


Figure 2 
Computing Circuits Utilizing the Property of a Vacuum Tube 


The A.C, impedance in the cathode circuit is great enough to 


cathode follower. When series 
use constant current devices 
potentials to the tube. In prin- 


ciple the CCD does not affect the voltage drops across R,C,C; 
and L or the currents through them. Pentode circuits, vacuum 
photo cells, and transistor circuits canbe used for this purpose. 
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Figure 3 

Typical R-C Feed-back Computer 
The A’s are high-gain inverting amplifiers. Similar circuits 
are possible for the other integro-differential equations. Suc- 
cessive differentiation by this method is to be avoided if pos- 
sible. High-frequency noise becomes emphasized and very iow 
noise levels are necessary in the power supplies and the trans- 
ducer. 


damped accelerometer to a step function in force, for example, 
would be a jump superimposed on a decaying sine wave. We can 
excite this response experimentally by applying a constant force 
through a string and cutting the string. In a series of repeated tests, 
the computer is adjusted until it converts this oscillatory signal into 
a step function with a very sharp rise time. 

More elaborate adjustment procedures are also possible and 
may be preferred wher the appropriate input functions are not prac- 
tical or when the required settings are numerous. We can then 
establish the relative values of theparameters of the transducer by 
means of the response to sine waves of various frequencies. This 
procedure has been developed for the analysis of servomechanisms 
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and other linear systems. There is also an easier approach for 
some kinds of transducers. The required data are available from 
the detail of the response to step functions and unit pulses. In these 
methods we derive numerical values descriptive of the transducer. 
The computer can, therefore, be designed accordingly or the cir- 
cuits can be adjusted exactly. 

There are two kinds of sensing elements that we have not dealt 
with in our examples of seismic transducers. Since these two are 
commonly used, it appears desirable to discuss the relationships 
which apply. These sensing units are the piezo-electric crystal and 
the inductive pickup. The inductive pickup produces a signal which 
is proportional to the velocity difference between the case of the 
transducer and its internal seismic mass. In the piezo-electric 
method, the displacement between the mass and the case causes a 
charge to be developed across the faces of acrystal. This crystal 
also acts as an electrical capacitor. The charge tends to leak off 
of this capacitor through accidental paths of conductivity and through 
the resistors of the circuits to which it is coupled. This leakage 
effect results in distortions for motions that are slowly varying. 

First, let us consider the inductive devices. The output voltage 
Eo is proportional to Y. Again we will assume that the constant of 
proportionality is 1. Consequently, the equations relating output to 
input are 


- eee - wu: K 
(7) X= -(E+ 7 E+ 7 ©), 
“ =F K 
(8) X = -(E + M E + 3 J Bae, and 
, u K 
(9) K = -(E + i [Bat + FF ff Bat at. 


The third derivative of the displacement (jerk) is included be- 
cause it can be approximated directly by the inductive pickup. On 
the other hand, the displacement itself cannot be approximated 
directly and is not listed. The conditions of the approximations are 
similar to that which holds for equations (2), (3) and (4). Equations 
(7), (8) and (9) are the relationships which are to be simulated in the 
analog computers when we attempt to reconstruct the original input 
to the inductive transducer. 

The pertinent property of the piezo-electric crystal is that the 
rate of formation of charge is proportional to the rate of distortion. 
If q represents the charge, and as before Y represents the displace- 
ment, then: 


(1¢ 


lea 
the 


ono &@ 
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(10) q = aY. 


When we take into consideration the formation of change, the 
leakage resistance (r) and the capacity (C) of the crystal, we have 
the system which is shown in abstract form in Figure 4. 


|-_—» 





al CURRENT 
GENERATOR 








Figure 4 
Equivalent Circuit Piezo-electric Transducer 


The generator produces a current which is proportional to the 
rate of change of deformation. The charge tends to collect on 
the capacitor, but is partially bled off through the leakage 
paths, R. 


According to Kirchoff's laws 


(11) q = -i +aY¥,, and 
Eo ° 
(12) CEo . R Fay. 


Once again we can let the coefficient of proportionality (a) be 1. We 
combine equation (12) with (2), (3) and (4) to determine the correct- 
ing equations for the piezo-electric transducer. The equation for 
the third derivative will also be included since it too may be approx- 
imated directly by the transducer for the proper combination of 
operating conditions and parameters. 


a KC | HO K 
(13) X -[ce+(4 + ar +e &+ sh | 


(14) & = [ ce + (4+2)8 + (AG. e+ am / Eat | 
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(15) Xx | cE + (4 +2 )z re + fsa +5 feaat | 


(16) x = -[ ce+(5 +2) fea + (% EC +H) fearat 
2 z Sf Edtdtdt | 


The computing circuits provide a method which in theory is ca- 
pable of extending the useful range of transducers as faraas we 
would like. In practice we have to recognize certain limitations 
First of all we rarely can determine an equation which is completely 
exact in representing the behavior of an actual transducer. What is 
more, the electrical analogs are not perfect either. Finally, we 
may have to compromise the entire effort by the use of a filter to 
lower the background of unwanted signals. Nevertheless, through 
the use of this technique, we have been able to extend our application 
of transducers to an &ppreciable extent. 
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Simplified Numerical Solution of a Catenary 


By Elihu Geer 
Professor of Civil Engineering 
University of Detroit, Detroit 


SYNOPSIS 


Methods for the solution of the several cases of a catenary are 
presented, with examples, the only tables needed being those of com- 
mon logarithms, and excellent accuracy being obtained with only one 
or two trials. The methods should be useful in solving problems 
arising in connection with tape measurements, cable conveyors, and 
electrical transmission lines. 


The following difficulties are met in using the usual formulas 
for numerical solution of a catenary: 

1. The solution involves cut and try, with the unknown both inside 
and outside a hyperbolic function or exponential. 

2. The tables of hyperbolic functions which are commonly avail- 
able are too coarse to yield accurate results. 

3. Initial trial values are not apparent, and convergence may be 
slow at first. 

There is offered below a system for the solution which has the 
following advantages: 

1. In certain cases the initial trial value may be of satisfactory 
accuracy. 

2. In many cases interpolation or extrapolation from a pair of 
trial values may yield a satisfactory answer. 

3. In very unfavorable cases only a second interpolation may be 
needed. 

4. No other tables than common or natural logarithms are used. 

The notation for the simple case is shown in Figure 1. 
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C is defined as ”, or Tp = C.q 


q = load per unit length of cable 
s = length of cable from low point to support 


The derivation of formulas (some of which are well known) 
follows. From the condition that the cable is flexible and that con- 
sequently the slope of the cable equals the ratio of vertical to hori- 
zontal components of stress, 


(yy) YG - 8, ag = Yi+(%) -ax = b Ver oc ax 


dx q.C S 
- < x/ . -x/C) a . x 
(2) whence s = > (e~C -e ) = C.sinh c 
— = ees | “ot 
Since , C sinh G? by integration: 
(3) y = C.cosh a 


(4) Squaring and subtracting, y? = s* + C?, ory= Vs*+C?. 


But y = f + C whence 


1 
(5) y? = (fC)? = s* + C’, or C = 5 « (8° - f°) 


Inverting s+ y=C (sinh > + cosh = =CeVC , 


2 2 
(6) In = = & orx=C.n 2-5 =f in 


Cc 2f s-f 
or x _ + (s-f - (2.302585) - hog (s+f) - log (s-] 





2f 


The elastic elongation of the cable between low point and point of 
support = 


_(Tds _ faqyds _ favs? +c ds 
48*)",e *J AE” AE wad 


(7) As= ae [sy + c.x | 


stre 


rh 
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in which s, y, C are after elongation. A small amount of cut and 
try may be needed in this computation, but initial s should yield 
practically the final A s. 
CASE 1. Given: s, f, q, A, E 
Required: a, Tp, Ta, As 
Procedure: Solve (5) for C 
Solve (6) for a 
Tp=4.C 
TA = q-y = q(f+C) 
Solve (7) forAs 


Example 1. A tape which has a length of 100.000 feet under 
stress, is suspended with a sag of 1.000 foot below the ends. 


q = .017 Ib. ft., 
A = .005 sq. in., 
E = 28x 10° psi. 


Required: Horizontal distance between end marks 
Pull at ends 





Elongation. 
Solution: C = (50+ 1)(50- 1) = 1249.5 ft. 
2x1 : 
; 50+1 _ 
a = 1249.5 In 50-1 - 49.987 ft. 


Horizontal distance between ends 99.974 ft. 


Ta = .017 (1 + 1249.5) = 21.26 lb. 
.017 
2x.005x28x10° 


- 0076 ft. for half tape, 





Elongation = ° |(50x1250.5) + (1249.5x49.987) 


.0152 ft. for entire tape. 


Evidently the tape should have been 100.011 ft. long in an unstressed 
condition to give a horizontal measurement of 100.000 feet with a 
sag of one foot. 


CASE 2. Given: a,f,q 
Required: s, Tp, Ta 
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Procedure: Solve for length of parabola of same span and 
sag, and consider approximately equal to that 
of catenary, 


Fa (t} -20(t)' ta (ty -.. 


3 5 7 
os 2 2 camer 9 
a Vi af 2 init V1 af 
(9) 2 *\a ‘s* a” *\a 
Using this approximate s, solve (6) for x. In this computation, take 
two values of s to the nearest whole foot in the vicinity of the ap- 
proximate value, and interpolate or extrapolate to determine the 


value of s which will yield x = a. In case extrapolation is needed, 
repeat the process with the improved value of s. 


Then solve for C in (5) 
Tp - q.C 
Ta = q.y = q(f+C) 


—_~ 
co 
— 
Dm 
nu? 


w 
te 


Example 2. 
Given: a = 90’, f = 60’, q = 0.9#At., A = .2647°", 
E = 30x10° psi 
Required: s, Tp» Ta, elongation 


=a ~ 90 \i+ (2x80)" , 90" 2x60 2x 
Solution: s= » et) * 4x60 ee Ve + (2220) | 





= 112.1 ft. 
for s = 112, x= lL (ii2° -60 *) In 172 . = 89.1605 
2x60 52 
| or 173 _ 
for s = 113, x= 5x60 (ii3° -60 ) In 53 = 90.3910 
.8395 _ 
Interpolating, s = 112 + 1.2305 * = 112.68 
C= _ (rix68" -60 ) = 75.81 ft. 


2x60 


Tp = 68.2%, Ta = 122.23*, 


CAS 


(10) 


R 


it 


ie 
i, 
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___ 0.9 
4s 2x. 2647x30x10° * [(112.68x135.81) + (75.81x90) 





= .00125' for half span. 
Example 3. Very deep catenary: 


Given: 200’ span, 200’ sag. 


= £0 Vs.( zxa00)" , (100r (100)? in| 25200 , VA. fax200)* | 
. 100) * 4x200 100 100 





= 232.34 ft. — 
- / 232 -200 432 
for s = 232 , x= = ~9x200 hb = 30° = 89.74 
a=<n3 —x2,2 
? _ 233 -200 433 _ 
for s = 233 , x = —~9x200 In 33 = 91.96 
100-89.74 
Improved s = 232 + 91.96-89.74 ~ 236.62 
for s = 236, x = 97.87 
a s = 236 + ae ie = 237.10 
for s = 237, x= 99.80) =} ‘ 
100-99.80 
for s = 238, x = 101.72, Improved s = 237 + 101.72-99.80 


= 237.10 


Even for such an extreme case, the third trial was unnecessary. 


CASE 3. Given: s,a,q 
Required: ft, Ty, Ta 


Procedure: From the first three terms of the expression 
for half of a parabola, 





(10) f= = V5- TIS - 905 
a 


Use this f as initial trial value in (6) 
Solve for C in (5) 
Tp =q.C, Ta = @.y = q(f+C) 


Example 4. A wire 250 feet long is suspended from two points at 
the same level and 220 feet apart. Compute the sag. 
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~ 128 : = 
Solution: f = —= V5 125 _ Cc 
16 -//115 - 90 iio = 61 ft. 
186x64 , 186 c 
x 
Try f= 61, x = 2x61 In 64 = 104.098 | . 
_ 176x74 176 | 
Try f= 51, x = 2x51 In 74° 110.630 Cc 
4 10x.630 | Cc 
Improved f = 51 + “6.532 51.96 
_ 177x773 177 | CASE 
Try f= 52, x = 2x52 n 73 110.039 Figui 
" 039 
Use f = 52 + 591 = 52.066 
, 177.066x72.934 177.066 : 
Try f = 52.066; x °x52.066 In 72.934 110.00; OK . 
CASE 4. Given: 6, & Fas 
Required: a, f 
} 
Procedure: y . Ta 
q G 
C =) y’ -s° ' R 
f=y-C r 
Solve (6) for a = x 
Example 5. A 100 steel tape weighing .02 lb./ft. is given a 15 | 
lb. pull. Compute the sag and the distance measured. 
Solution: y 3 750 ft. 
c = ¥ (750)* -(50)? = 748.331 ft. P 
f = 750-748.331 = 1.669 ft. . 
Ne 
‘ ~. tn StY - 800 . 
a=x =C-ln C 748.331 In 748.331 49.963 ft. , 
Horizontal distance = 99.926 ft., neglecting stretch. I 
.02x3.4 r 


As 


. 2x.02x28x10° | (50x750) + (748.33 1x49.026) | .005 ft. for 


half tape 
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Corrected s = 50.005 ft. 











} 
: Corrected y = 750x a. 00 750.075 ft. 
Corrected C = |/(750.075)* -(50.005)? = 748.406 
800.080 
Corrected ; } 748. es 
a x 48.406 In 748.406 49.968 ft. 
Corrected Horizontal distance = 99.936 ft. 
CASE 5. The points of support are at different levels, as shown in 
Figure 2. 
‘ an ae Si 
nN " Sh Ps 4 
e:_™ , 
é re nen ge ; 
Figure 2 
Given: (a+b), fa, fg, q 
' Required: a, b, Sa, Sp, Tp, Ta, TB 
’ ~ _atb 
Procedure: From the properties of the parabola, a lei 
‘| B 
fs 
15 . 


Compute sa from (8) or (9). 
Choose a whole number of feet as a trial value for Sa, and com- 
SA -fa _ Sh - fB _¢ 


pute the corresponding value of sg in 
2fa 2fp 





(sa+fa) (sp+f 
Solve for span in (a+b) = C. In a Sp + fp) 
(s, -f,) (Sp - fg) 
Take another value of s, and repeat the process. 


Interpolate or extrapolate to obtain an improved value for sq and 
repeat until the derived value of a+b satisfies the given condition. 
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Figure 3. 


Example 6. The problem is shown in Figure 3. 


~ 600 _ _ 
“a vA 380.4 ft. 
150 





_ 2 i590 _\* 2 150 ) 4 pee ‘ 
sa = 380.4 + © 380.4 (=e) - £ 380.4 (= + > 380.4 ) =... 


416.7 ft. 
(416)° -(150)* sR -(50)* 





Try sa = 416, , Sp = 229.5, C = 501.85 





2x150 2x50 
? 566 279.5 
e +} ° — . 
a+b 501.85 - In 266° 179.5 601.18 


(415)* -(150)* __s*p -(50)* 





, Sp 7 228.93, C = 499.08 





2x 150 ~ 2x50 
? . 1p, DOD. 278.93 os 
a+b = 499.08 - In 265 178.93 999.43 
57 
Improved value of sa = 415 + 75 415.32 


Try sq = 415.32 ft., sp = 229.12, C = 499.97 ft. 


565.32 279.12 
265.32 179.12 





a+b 2 499.97 - In 599.984 


Sa = 415.33, sp = 229.128 C = 499.997 ft. 


? 3 279.128 
a+b = 499.997 - In 565.33 , 279.128 


265.33 179.128 600 ft.; OK . 
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CASE 6. Given: (S,+Sp) = M, (f,-fp) = K, (a+b) = L 
Required: a, b, fy, Sa, Sp, fp. 
Procedure: Assume a value for sa. 
Compute the corresponding value of f, and C in 


— f, _ (M-s,)* -(f4-K)’ 
2f, 2(f, - K) 





Try these values of sa, fa, and C in 


a+b = C. In (sa+f,) (M-s q+fa-K) 


(s avfa) (M-s y-f a+ K) 





Assume another value for s,s and repeat the above. 


Interpolate or extrapolate to an improved value for sa and repeat 


the process until a+b satisfies the given condition. 
Compute remaining unknowns as before. 
CASE 7. See Figures 4 and 5. 
Given: q = weight per unit length of cable 
W = concentrated load 
K = difference in altitudes of supports 
L = Total span 


" s', = length of cable from concentrated load to high 
support 
S'p = lengthof cable from concentrated load to low 
support 
Required: Location of low point, tensions in cable. 


Procedure: Due to equilibrium of horizontal forces, 





(11) c - Safa - Sp fB - si - ff - 82 -f 
2f 2fp Qf; of, 
Since (f, - fi) - (f, f.)=K, 
(12) Vs. +C? - Vs, +C* - Vsi + C? + Vs? + C? = K 
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Due to equilibrium of vertical forces, 


(13) Si + S2 ~ 
q 


W 
J=s'’a+8'RB+— 
q 


(14) sa + Sp = (Sk +581) +(S’p +s 


tS 


Assume one of the four lengths and compute the others. A reason- 


W 
able starting point is s; = S2 2q° In that case, 


s’, +K* -s},)’ a 
, S -c* 
c=} 2K he 


Solve for C by cut and try in (12). 








Then f, = Ys, +C° -C, fp = Vs%+C? -C 
f, = Vsi+C? -c, f, = Vs? +C? -c 
a C. ln SA +fa b =C.ln sB+IB 
SA- fa Sp -fB 
x, = C. ln S: +h te =C. in? fa 
4° f, S2- f, 
Then test to see if a-x; + b-x, = L. 


If not, try another s, repeat, and interpolate. 

The work in this test may be shortened somewhat by using: 

(15) L C. In (s, +fa)(Sq +f p)(S1 -f1)(S2 -f2 ) 
(Sa -f \(sp-fp)(si +f, )(S2 +f, ) 





instead of computing a, b, xX,, X2 


2 individually. 

Example 7. The problem is given in Figure 6. 

600 
10 


Solution: $1+Se 60. Assume s; = Sz = 30’ 


Then s, = 209 + 30 = 239’, sp = 105’ + 30’ = 135’ 











ELIHU GEER 


f; = f, = 1.007’ Xi = Xe = 29.968’ 
f, = 59.974’ a = 228.837’ 
fg = 19.974" b = 133.123’ 


a+b-x:-x2 = 302.024’. 


< L=300 - 







q=-10 Ib/Fft 






Ss 


y W=600 Ib 


Figure 6. 


A fair approximation for the location of the low point now can be 
found by assuming that the chord lengths from low point to points of 
support remain constant. 

The horizontal location thus computed is only .04 foot in error, 
while the vertical location is approximately 0.9 ft. too low. There 
is no obvious choice for the next trial values for lengths. Trying 
SA= 235’, sp= 139’, s,; = 26’, and sz = 34’, it is found that 
C = 412.93’, fa = 62.187’, fp = 22.767’, f: = .817, f, = 1.397’, and 
that a+b-x,-x, = 300.457’. Extrapolating, a value of slightly less 
than 234' for sa is indicated. Trying sa = 234’, sp = 140’, s, = 25’, 
Ss, = 35’, it is found that C = 404.560’, fa = 62.799’, fp = 23.539’, 
f, = 0.772’, f2 = 1.511’, a = 222.597’, b = 137.346’, x: 24.984’, 
X. = 34.957’, which satisfy the conditions of the problem quite well. 
The 600 pound load is 102.39 feet from the left support measured 
horizontally, and 22.03 feet measured vertically. The tensions are: 


Ta = 10(62.799+404.560) = 4673.59 Ib. 
Tp = 10(23.539+404.560) = 4280.99 Ib. 
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The Mathematical Background of Spiroid Gears 


By Oliver Saari, Development Engineer 


illinois Tool Works, Chicago 


1. THE GENERAL FIRST-ORDER GEAR TOOTH ACTION LAW 


Spiroid gears are more easily understood if some basic concepts 
regarding gear tooth action are first clarified. In fact, it was the 
formulation of these concepts which led to the development of the 
new gears. 

Most gears rotate on fixed axes, at a constant angular velocity 
ratio. The axes can be parallel, intersecting, or skew, and many 
types of gears have been made for operation in each of these relative 
positions. Conjugate action of all of these gears can be shown to be 
in conformity with a perfectly general first-order gear tooth action, 
which will be here briefly outlined. 

In general, active gear tooth faces are continuous surfaces having 
continuous first derivatives, i.e., they are smooth surfaces having 
no gaps or folds. They have instantaneous tangential contact at one 
or more points, which will be called “conjugate points.” If the gear 
teeth have line contact, which is most often the case, there are an 
infinite number of these conjugate contact points for any instantane- 
ous position of the gears, and the aggregate of these points makes 
up the contact line. The requirement that gears transmit motion 
between two fixed axes at a constant angular velocity ratio imposes 
certain restrictions on the tooth surfaces at every one of their 
conjugate points. 

It is convenient to think of the contact at a conjugate point as 
being two infinitesimal particles, lying instantaneously in the com- 
mon tangent plane. In general these particles have a relative motion, 
which results from the rotation of the gears on their axes. The 
fact that one gear tooth is imparting motion to the other implies that 
the direction of relative motion of the contacting particles must lie 
in the common tangent plane of the tooth surfaces, as illustrated 
in Figure 1. 

The absolute velocities of the contacting particles are V,-Vo. 
Their relative velocity is the difference, Vi-V2.. From the laws of 
vector subtraction, this is known to be a vector parallel to the line 
ac which joins the ends of Vi and V.. From study of Figure 1 can 
be inferred the fact that if the relative velocity lies in the tangent 
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Figure 1 


Basic Velocity Relationships 


plane, then the normal components of the absolute velocities are 
equal, and conversely. This is true in any case where one surface 
drives another through a tangential contact point. This, and the fact 
that gears must operate on fixed axes at constant angular velocity 
ratio, are the conditions which determine the tooth-action law. This 
law is called “first-order” because it deals only with the tangency 
and common tangent plane of the tooth surfaces and says nothing 
about higher-order properties such as surface curvatures. A fa- 
miliar first-order tooth action law is the two-dimensional law for 
transverse sections of parallel-axis gears: the common normal 
to mating tooth curves at every conjugate point passes through a 
fixed point, called the “pitch point.” It is not surprising to find that 
the general first-order tooth-action law includes this as a special 
case. 

The most general disposition of axes in a fixed frame of refer- 
ence is illustrated in Figure 2. 

For consistency in nomenclature, let us call the gear axes #1 
and #2, with axis #1 always so placed that it coincides with the x- 
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Figure 2 


The Most General Disposition of Axes 


axis of the Cartesian coordinate system. The angle between the 
axes is @ and the shortest distance between them is C. The axes are 
parallel when a = 0. They intersect, as in bevel gearing, when 
C #0 but a#0. They are skew, as in worm, hypoid, and spiroid 
gears, when C # 0 anda#0. 

This general case will be returned to later. For simplicity in 
demonstrating principles, let us first consider the right-angle skew- 
axis case, for which C is arbitrary and a = 90°, as shown in Fig- 
ure 3. Moreover, let us restrict the analysis to a left-hand system, 
in which the angular velocities of the gears, W, and Wz, are directed 
along the axes as shown. Let K be the fixed angular speed ratio of 
the gears, and p (x,y,s) be any conjugate point. At every conjugate 
point p, then, are two conjugate particles, each fixed to and rotating 
with its associated gear. The location vectors of these particles are 
f; and Ff, and their absolute velocities are v; and V2. 


Let 
N, = number of teeth in gear #1 
N2 number of teeth in gear #2 
K angular speed ratio = N2/N, = |W,l /|We 


Let i, j. and k be the usual unit vectors associated with the 
X,y,Z axes. We may assume any convenient value for the angular 
speed of one gear, as long as the speed of the other is then made 
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Figure 3 
Right-Angle Skew-Axis Gear System 


consistent with the speed ratio K. The following values are con- 
sistent: 


Wi -Ki 
— a 
W2=) 


The location vectors of the conjugate particles can be seen from 
Figure 3 to be 
Ti xi + yj + zk 
Te = fi - Ck 


Absolute velocities are given by cross-multiplication according 
to the rules of vector analysis, 


Vi = Wi Xf = (Kz)j - (Ky)k 
V2 We x T2 = -(C-z)i - xk 


Let 


relative velocity = V; - V; 
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(1) V=(C-z)i+ (Kz)j + (x-Ky)k 


The components of the relative velocity vector are functions of 
the coordinates. Obviously its direction at any point is not affected 
by speeding up or slowing down the rotations of the gears, as long 
as the angular speed ratio is maintained. Equation (1) defines a 
vector field of infinite extent, since obviously every point of space 
referredto the x, y,z coordinates has some value of relative velocity 
associated with it. Such a vector field is more easily visualized if 
we can discover its “field lines,” or the family of space curves 
which are everywhere co-directional with the vectors of the field. 
The differential equation of the field lines is 


dx : dy : dz=(C-z): (Kz) : (x-Ky) 


This can be more easily solved if coordinates are transformed 
as follows: 


a 

x V1+K? (Kx -y’) 

y Vi+K" (x +Ky ) 
ee = 
+ LK 


This amounts to a rotation and translation of the axes to a new 
position shown in Figure 4. 

The translation and rotation of the coordinate axes are deter- 
mined from the speed ratio K, thus 


tan B =— 


nin 


Sins 


h +K’ 


— 


Let i’, j’, and k’ be unit vectors corresponding to the new co- 
ordinate axes. The transformed equation of the relative velocity 
vector field becomes 


Q) V- Vink Cea z'}'- ye| 


and the differential equation of the field lines becomes 











OLIVER SAARI 


1 AXIS OF GEAR #2 








¥ AXIS OF GEAR #1 


Figure 4 


Transformation of Coordinates 


o-« 0 « —_— -« , 
dx’: dy’: dz’ = Dor: 2’: -y 


This system of simultaneous differential equations can be easily 
solved, resulting in the equations of two surfaces whose intersecting 
curves are the field lines, 

2 


y” + z” = an arbitrary constant = P* 


» (By oie’) _ : 
x’- OK? sin +) = another arbitrary constant. 


The first of these equations represents a family of cylinders whose 
axes coincide with the x’-axis, and which have arbitrary radii Pp, The 
second equation represents a family of left-handed radial-helicoidal 
surfaces of equal lead, axes on the x’ axis, and arbitrary phase re- 
lationship. The lead LR is given by 


Lr = 2286 


1+K 


Simultaneously the two equations represent all left-handed helical 
space curves of the above lead, whose axes coincide with the x’- 
axis. These helices are the easily-visualized field lines of the 
relative velocity vector field. The x’-axis (Figure 4) will be called 
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the relative axis and the field lines the relative helices. For the 
right-angle gears, giving the center distance and speed ratio deter- 
mines the position of the relative axis and the constant relative lead. 
For a given center distance and ratio, it is evident that through 
every point in the space passes one and only one relative helix. The 
relative axis is itself one of the relative helices, corresponding to 
p= 0. 

These methods can be extended to the more complex case of 
general disposition of the axes, where a, C, and K are all arbitrary. 
In terms of these quantities the position of the relative axis and the 
value of the relative iead are 


sina 
7 ase K+cos a 
C(1+K cosa 
(4) a- 1+K* + 2Kcos a 
27KC sina 
= . 6 6 a 
(5) LR 1+K" + 2Kcos a 
These are illustrated in Figure 5. 
The general first-order tooth-action law can now be stated as 
follows: 

A necessary condition for conjugate action at any point is that the 
common tangent plane must also be tangent to the relative helix at 
the point. 

That this is a perfectly general law, reducing in special cases to 
the well-known laws of conjugate action, can be readily demon- 
strated. For example, in parallel-axis gearing, since a= 0, 


B =0 

Co 
h = iK 
Lr = 0 


Here the relative axis is a line parallel to the axes of the gears and 
co-planar with them. It is the tangent line of the pitch cylinders as 
conventionally defined. The relative helices have zero lead, being 
annular circles concentric with the relative axis; therefore in any 
transverse plane they appear as concentric circles with center at the 
pitch point. The general law is here equivalent to stating that the 
trace of the common tangent plane must be tangent to the relative 
circle passing through the conjugate point. But, as illustrated in 
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Figure 6, this is exactly the same as saying that the common normal the } 
of the contacting tooth profiles must pass through the pitch point. can 
For intersecting-axis gears, C = 0 anda+#0O. Substituting these point 
values in equations (3), (4), and (5) gives C 
| easil 
tan Bp = Sno cons 
K+ cos a helix 
h =0 h 
oluti 
abou 
Here the relative axis is the contacting element of the well-known 
bevel gear pitch cones, and the relative helices are again annular 
circles about this line. Here the general law is seen to be equivalent 
to the statement that all common tooth normals at conjugate points q 
must pass through the pitch line. surf: 
In these simpler cases, no apparent advantage results from re- mal : 


stating the tooth-action law in terms of relative helices instead of | 
having tooth normals intersect the pitch lines. It is in its applica- 
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Figure 6 


General Law as it Applies to Parallel Axes 


tion to skew-axis gearing that the general law gives insight not an- 

ticipated by the conventional ways of stating the law. In skew-axis 

gearing the common tooth normals do not necessarily pass through 
‘mal the relative axis. At any conjugate point, the common tooth normal 

can lie anywhere in the plane normal to the relative helix at that 
ese point. 

Certain interesting facts about the relative helices can be fairly 
easily proven. The amount of “sliding” action between gear teeth is 
constant alone any relative helix and is always directed along the 
helix. It is a minimum along the relative axis itself. 

Moreover, any relative helix generates tangent surfaces of rev- 
olution about the axes of the gears. For example, if a given relative 
helix is used as a generatrix to describe a surface of revolution 

} about axis #1, a surface normal vector is given by the formula 


nown 

vular Vv = Vi x (V; - V2) = = v; x V2 

alent ni 

oints If the same relative helix is used as a generatrix to describe a 
surface of revolution about axis #2, the corresponding surface nor- 

1 re- mal vector is given by 

ad of 


lica- 
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—_ == al ao => —_ 
ae = V.x (V; - V2) > -V, xV2 
But these surface normal vectors are identical. Therefore the 
surfaces of revolution are tangent to each other at every point of the 
generatrix. The surfaces have constant sliding velocity at every 
point of contact, this velocity being directed along the generatrix. 
One such pair of tangent surfaces is the pair of hyperboloids of rey- 
olution generated by the relative axis itself, the well-known “pitch 
surfaces” of skew-axis gearing. However, surfaces of revolution 
generated by any of the other relative helices have the same prop- 
erties and might just as well be called “pitch surfaces.” 


2. Application of the General Tooth Action Law to Skew-Axis 
Gear Design 


Most skew-axis gears are produced by making the teeth of one 
member of some arbitrary form and generating the teeth of the other 
member by some cutting process which duplicates the rolling action 
of the gears. Conventional worm gearing, formate hypoid gearing, 
and spiroid gearing are of this type. In each case, the tooth form of 





— 


- — 


the primary member is first chosen as one which can be made by | 


some simple form-cutting process. A cutting tool is made whose 
form represents as nearly as possible the teeth of the primary gear 
member. This tool is then rotated with the work-blank of the con- 
jugate gear member in a manner which duplicates the relative motion 
of the gears in action. 

This method results in useful gears only if the tooth-form of the 
primary member is chosen with some care: it must be such that 
conjugate action is possible over a large portion of the available 
tooth surface. An analysis of possible contact lines in this type of 
gearing can be made by applying the general tooth action law as fol- 
lows: 


Let 
(6) F(x, y,z) = 0 


represent the pre-determined tooth form of gear #1 in some in- 
stantaneous angular position. Direction numbers of the tooth normal 
at any point are the partial derivatives F,, Fy, and Fz. For con- 
jugate action, the tooth normal must also be normal to the relative 
velocity vector. Therefore 


~ 


(7) 
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(7) ¢ (x,y,z) = 0= VF, + VyFL + VZF, 


yy 
where Vx, Vy, and Vz are the components of the relative velocity 
vector V. 

If the surfaces defined by equations (6) and (7) intersect, their 
intersection is the instantaneous contact line. If they do not inter- 
sect, conjugate action is not possible. Interesting general conclu- 
sions can be drawn from an analysis of this sort. 

For example if the arbitrary tooth form of gear #1 is a plane 
surface, then F!>,y,z) is of no higher than the first degree in any of 
the variables. ’<s partial derivatives are therefore constants. Since 
the components of the relative velocity vector are also of first de- 
gree, ¢ (x,y,z) is also a plane surface. The contact line is the in- 
tersection of two planes, which is always a straight line. This is a 
perfectly general result which can be proven in many more laborious 
ways. 

The contact analysis is not always as easy as this, particularly 
when the function representing the primary tooth surface is in para- 
metric form. However, the zone of possible conjugate action can 
always be found explicitly point by point. 

Given a form-cutting method by which we wish to produce the 
teeth of the primary gear member, there will always be some arbi- 
trary parameters. In spiroid worms, which are cut with a taper- 
helical motion of constant lead, we can vary the lead, taper angle, 
position, and the pressure angles of the teeth. The tooth action law 
helps us choose these parameters in a way which will insure con- 
jugate action. 

Referring now to the specific right-angle skew-axis case, Fig- 
ure 3, let us say that conjugate action is desired at some arbitrary 
point (called the “pitch point”) which is designated by giving numer- 
ical values to its x,y,z coordinates. If the center distance C and 
speed ratio K are given, then the relative helices are determined, 
and one of them passes through the pitch point p. Its direction in 
space is the same as that of the relative velocity vector given by 
equation (1). 

It is convenient to relate the relative velocity vector to a new 
Cartesian coordinate system whose origin is at the arbitrary point 
p, as shown in Figure 7. 

Point p is the arbitrary point, and r is the radial distance from 
axis #1 to this point. The new coordinate system is designated by 
three mutually perpendicular vectors 1, 4, and b, originating from 
p and defined as follows: 


i is the x-unit vector, and is therefore parallel to axis #1. 
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Figure 7 


Transformation of Coordinates 


‘ 


a is a unit vector in the direction of the radial line through | 
p, and is therefore perpendicular to axis #1. 


* b is a unit vector at right angles to both i and a, defined by 

b =4xi, and is therefore at right angles to both axis #1 and the} 

radial line r. 

These three new unit coordinate vectors can be expressed in | 
terms ofthe x, y, z coordinates of the point p as follows: 


_—_ 


1=1 


© inser = 
a = — (yj + zk) 


- 


—_ l : a 
b == (2) ~ yk) 


Let Vx, Va, and Vj), be the components of the relative velocity 
vector in these coordinate directions. Then 
al 


Vx =V. sC-g 
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———— XZ 
Va =V.a a 


Vp = V.b = Kr - = 


The relative velocity vector at p is tangent to the surface of rev- 
olution (or pitch surface) generated by the relative helix through p. 
Let tT be the taper angle of this surface at the point p, as defined in 
Figure 8. This instantaneous taper angle of the pitch surface at p is 
given by 
V 
(8) tanT —— , — 

Vx r(C-z) 
The relative velocity vector is tangent to a cone of this taper angle, 
whose axis is on the x-axis, and which contains p. Hence it is co- 
directional with some helix of constant lead drawn in the cone. This 
helix must have a lead L which is related to the axial and tangential 
components of the relative velocity vector as follows, 


Vx 2n(C-z) 
(9) L = (27r) % hes 
K -~3 
x 
= z 
| 
TANGENT 
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Pitch Surface and Tangent Cone 
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Equations (8) and (9) are fundamental formulas in spiroid gear- 
ing. A spiroid worm which has a taper-helical element of this taper 
angle and lead in a conical reference surface containing p will have 
possible conjugate action at p regardless of other tooth form para- 
meters. In particular, if both sides of the worm thread are gen- 
erated along this tapered helix, they will have conjugate action at p. 
The amount of conjugate action on the remainder of the thread sur- 
faces depends on the pressure angles and on the location of the pitch 
p. In general it is better to have the pitch point on the outside of the 
worm blank as this minimizes the amount of inactive fillet generated 
in the mating gear. 

It is interesting to note that conventional cylindrical worms are 
but a special case of the gearing defined by these formulas. Here 
the pitch point is on the common perpendicular between the axes of 
the gears (the z-axis). Thus x= 0, y= 0. Substituting in equations 
(8) and (9) gives r= 0 and L = 27(C-z)/K. But this means the pitch 
surface is cylindrical and the lead is exactly that of a conventional 
worm whose pitch radius is z. Furthermore, if x =0, y #0, and 
z= constant, we get the same taper angle and lead from the for- 
mulas. This means that pitch points chosen on a line in the y-z 
plane, parallel to the y-axis, all give the same worm thread para- 
meters. Therefore the zones of contact for both sides of the worm 
threads intersect along this line regardless of pressure angles. But 
this is simply the well-known pitch line of conventional worm gear- 
ing. 

Spiroid gears should be thought of as a generalized form of worm 
gearing. They apply the foregoing principles to the purpose of ex- 
tending and improving the nature of the tooth contact lines. 
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Design of Orifices for Hydraulic Shock Absorbers 


By T. E. Lohr 


General Motors Corporation 
Detroit 


The modern day automobile, as it increases in size, weight, and 
power, is becoming what would be classified ten years ago as a full- 
fledged race car. This type of vehicle for our present day highways 
imposes serious problems for the chassis engineer. Among the 
most pressing are: (1) high speed side sway control, (2) longitudinal 
“dipping” control, and (3) isolation of the body from road bumps and 
noises. A great deal of development has been done in suspension 
systems such as the “air-ride,” the full and partial torsion ride, 
and the torsion air ride, to attempt to improve roadability. How- 
ever, it is apparent that regardless of which type of suspension sys- 
tem is used, the success or failure of that system depends either 
partly or wholly on the shock absorber selected. In the infancy of 
the automobile, shock absorbers were designed by trial and error 
until a suitable unit was chosen for the vehicle. Today, this field 
has become an exact science. 

Until recently, when a large manufacturer announced the use of 
anew aircrait type shock absorber known as the “oriflow,” we were 
still applying old concepts of design. This new method of ride con- 
trol opened the eyes of the industry and has led to even greater 
improvements in shock absorption. This article will deal with two 
basic types of orifices for future consideration in ride control 
developments. 


Symbols used in the computations: 


= Total orifice area, in inches’ 

= Piston area, in inches’ 

= Orifice discharge coefficient 

= Diameter of cylinder, in inches 

= Piston diameter, in inches 

Piston diameter at start of travel, in inches 

= Kinetic energy absorbed by fluid, in foot-pounds 
= Acceleration, in foot-pounds per second’ 

= Force on piston, in pounds 

= Acceleration of gravity = 32.2 feet per second’ 
= Pressure, in pounds per inch? 
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S = Piston travel from start, in inches 

Smax = Full piston stroke, in inches 

V = Velocity of piston, in feet per second 

V = Velocity of fluid, in feet per second 

W = Weight of body to be decelerated, in pounds 

y = Density of fluid in cylinder, in pound per foot* 

(hydraulic fluid = 69.4 lbs/ft*) 

M Mass W/g 

Pa eo: ee lan 
Ww TYPE 
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Actually, there are five types of orifices although only two are 


in widespread use. 


It is these two that we are concerned with. 





EQUATIONS: 

a) w = rca? dS /3 456 a? 

(2) P= re AY 

(3) Type l Log, = : VAS. 

(4) Type 2 WLoge _ cry mm +) |- ae . 
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To illustrate the application of the equations, a typical shock ab- 
sorber problem will be considered. The various orifice sizes will 
be calculated and the performance will be analyzed in terms of 
cylinder pressure and piston velocity versus piston stroke. 


Problem: 

A 1000 pound weight traveling at a rate of 100 feet per second is 
to be brought to rest in a distance of three feet. Space limitations 
require that the piston diameter be approximately six inches. 


TYPE 1 ORIFICE 


Assume: 
a = 3 inches* — Substituting in equation 3, we have: 


Loge (100/v) = 76 x 1.2? x (0.7854 x 6*)*s/(3456 x 3” x 1000), 


from which it is possible to solve for values of v for corresponding 
values of s. 


By using equation 2, also, the following tape has been prepared: 
s 0 2 4.5 9 18 27 36 
Vv 100 85.2 70 49 24 11.7 5.75 
P 10500 8260 5600 2520 604 153 334.7 


TYPE 2 ORIFICE 


Assume total orifice areas to vary between 8 inches’® and 0.25 
inches’ , divided among three orifices. The piston diameter is six 
inches, and the depth of slot is 0.75 inches. 

The central angle subtended by each orifice will then vary be- 
tween 60.4 and 1.9 degrees, for which the corresponding chords will 
be 3.4 inches and 0.111 inches. This will give K a value of 0.0915. 
By use of equations (4) and (2), the following table is computed: Note 
that in equation (4): 


T = Depth of slot = 0.75 inch 

N = Number of slots 3 

¢@ = Central angle subtended by slot radians 
K Aue where: 


Xo = Maximum width of slot 
X = Width at given station 
S = Piston stroke from initial position. 
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S 0 3 9 12 18 22.5 27 30 8636 
xX 3.4 3.13 2.59 2.30 1.75 1.34 0.93 0.66 0.11 
@ 60.4 55.2 44.4 40 30 25.3 15.82 11.2 1.89 
a 8.00 7.31 5.88 5.30 3.97 3.03 2.09 1.49 0.25 
Vv 100 93.2 85.7 78 64.6 50.7 24.9 12.1 0 
P 1477 1530 2000 2040 2570 2713 1365 645 0 


EFFECT OF FRICTION AND SPRING FORCES 


In the foregoing discussion, piston velocity is assumed to be 
affected only by orifice size. This is the ideal case assuming that 
friction is negligible; however, any or all of the following forces 
may exist: (i) mechanical friction, (2) restoring forces exerted by 
springs for returning the piston to zero position, (3) packing friction, 
(4) viscous friction, (5) mechanism weights when arranged in a ver- 
tical direction, and (6) hydraulic fluid elasticity. Viscous friction 
and fluid elasticity will be taken up later. The method of including 
constant and variable resistance will be considered here: 


Kinetic energy to be dissipated = Wv*/2g. 
Frictional energy (mechanical and packing) = K: s/12 
where K, = frictional load. 
Spring energy = {(So+S) K, ds /12 where 
So = pre-compression, in inches 
K, = spring constant, in inch/pounds. 


For a constant decleration orifice, merely subtract the frictional 
energy from the kinetic energy and solve for the necessary orifice 
area “a”. For all other cases, the procedure is as follows: 

The energy equation is 


5 2.243 . de * (So+S)K2 de 
wove - v*) 26 =f TOV Ads . ( oo ior oe 


3456 ga t” 


0 “0 0 


12 


Differentiate the energy equation above: 
-Wyd,/g = Ksv’A® d,/a® + (Ki+K2So) dg/12 + K,Sds/12 where 
Ks = YC,/3456 ¢g. 

Divide thru by dg: 


- Wy dy 
gd 





= K3v’A*/a” + (K, + K,S)/12 where K, = K; + Kz so let 
Ss 


v=Vy then = & , KaAay » Sarit : 


2g ° d, a 12 
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This equation is a linear, first order differential equation, the solu- 
tion of which is 


y=v’ =e-] K;2gA°d,/(Wa’) { - fe fks2ga° a, /wa? 2g(K, + K2S) 
d,A2w+e 


from which it can be seen that for S = 0, c = v’0. It may be noted 
that the value of Js 2gA; /(W,’ )\ds may be computed directly by 


multiplying the right side of (3), (4), and (5) by 2 W. 

In order to obtain an idea of the comparative effect of frictional 
forces, their effect on piston speed will be considered. 

Assume that there is a constant mechanical friction force of 250 
pounds due to a coefficient of friction of 0.25, also, that there is a 
constant packing force of 960 pounds. This value is determined by 
multiplying the highest pressure to be encountered by 0.2. There is 
a spring retarding force which is from 1000 pounds at the beginning 
of the stroke to 2000 pounds at the end of the stroke. Let the prob- 
lem be to determine the velocity of the piston at 18 inches or 50 
per cent of the stroke. By equation (6): 





Ss 

3 

(6) { Ks 2gA d,= aw Loge (v/vo) =-0.391. Also 
0 W Ww 


K, = Ki + Kz So = 250+ 960+ 1000 = 2210, 


substituting in equation (6) and solving step by step integration, 
v = 81 feet/second compared to 82.2 feet/second calculated without 
including frictional effects. In this case, the effect of friction is 
negligible, but it cannot be assumed that the difference will be neg- 
ligible in every case. Usually the friction effect may be neglected 
in making the first approximation for orifice size. 


EFFECT OF VISCOSITY 


Frictional losses due to fluid viscosity are small in shock-ab- 
sorber design. In the example chosen for this article, assume that 
the piston has traveled 2 inches. The coefficient of viscosity, 
# = 0.00003 pounds per second/foot.* Assume an orifice of type 2 
is used: 


Fy, = ua v/y = wa AyAay) where a = circumferential 
area of the piston 











150 


3+ 10° x 67x 2x 0.7854 x 6? iv 


then Fy 





144 x 0.75/12 ” 


F,, for a velocity of 100 feet/sec will then be approximately 0.5 
pounds, which is negligible compared to a hydraulic force averaging 
about 40,000 pounds. 


EFFECT OF FLUID ELASTICITY 


The average volume modulus of elasticity of water at an operat- 
ing pressure of about 4000 p.s.i.is 330,000 p.s.i. Also Ey = volume 
modulus of elasticity in p.s.i.; P = Pressure, in P.S.L; V = Volume, 
in ft.* 

If maximum pressure developed using type 1 orifice is 10,500 
P.S.I. and the volume of the fluid is 0.59 feet®, then V = 0.59 x 
10,500/330,000 = 0.0188 cubic foot compression. This amounts to a 
linear displacement (using 0.5 foot diameter) of 0.0957 feet. The 
work absorbed by this compression is 0.5 x 0.0957 x 10500 = 503 
foot/pounds, which amounts to 3 per cent of the original kinetic 
energy. However, when it is considered that the pressure is re- 
duced to zero by the time that the piston reaches the end of its 
travel, it may be seen that the energy of elasticity has been reduced 
considerably. 





Continuing: 
- ; 144p C*Vv’ _ C,Vv*Y 
Bernoulli’s Theorem: - y * Qe or * 388 2° 
2 2 
hence the force on the piston= F = Ay = ae while the work done 


by the piston during its stroke = E f Fd,/i2. 
0 


This equals the change in kinetic energy of the mass to be stopped. 
Hence: 


W(vo - v*) | z ¥C*V* Ads 
2g 7 0 3456 g 


By continuing equation - V = Ay/a, therefore, where vo = velocity of 
weight (W) at start of stroke 


8 
Wve - v*) _ ¥C* Vv" A*d, 
2g 0 3456 ga* 


Differentiating and dividing three by v’ 








also 


Inte 


w i 





of 








r described earlier in the article. 
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dv  YC*A°d, 


‘Vv -3456a" 


also, we have from above P = ¥C*A’vy?/288 ga’. 


, 


TYPE 2 ORIFICE 


The usual method of solving for velocity is to divide the stroke 
into a sufficient number of small segments and compute the pres- 
P = 2c’ A*y’ 

288 ga? ga? and assuming that 
Vv = Vo is constant over the whole segment. The deceleration then is 
computed using P = wf/gA. The velocity for the next segment is 
found, remembering that v* = v6 - 2fs/12. This process is repeated 
until the velocity is found for every segment. The disadvantages of 
this method are as follows; (1) The fewer the segments taken, the 
more the point of maximum pressure is moved toward the beginning 
of the stroke. By integration, the actual point of maximum pressure 
is found to be nearer the end of the stroke. (2) Any error at the be- 
ginning of the calculations will necessitate starting over at the 
point of error; moreover, the error is cumulative. (3) And many 
times, after all of the work of computation has been carried out, it 
will be found that the entire stroke has not been used efficiently, in 
that v becomes negligible too soon along the stroke (thus introducing 
higher pressures than are necessary) or that v is too high at the end 
of travel, causing shock absorption to be inadequate. 

If the width of the slot is described by the following expressions: 
(Xo - x)/S = K where K/2 also Tan @/2) and sin (9/2) = xd, = 
(Xo - Ks)/de and K, = Xo - dy sin (6/2) and ds = -d, cos (9/2)d¢ 
/(2K). Also, let 


sure over the first segment, using 


n = number of slots 

t = thickness of slot, in inches 

¢ = angles of slot, radius 
Then: a = de ¢tn/2. 


wdv 2c*A*dco cos (9/2)d¢ 
v  6912K (d,/2)* t?n’¢’ 





Substituting: 


Integrating we have: 


yc? A? | -cos(¢/2) (¢/2)° —(/2)° +4 








7 
WG. y, 13824 d,K(tn/2)” | (2) rss 5L5 


The remainder of the calculations would then be similar to those 
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